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I.  INTRODUCTION 


As  recently  as  1960,  very  little  was  known  about  the  basic  mechanisms  of 
fan  noise  generation  and  transmission.  In  the  absence  of  well  validated 
causal  relationships  between  parameters  and  noise  characteristics,  researchers 
in  the  field  first  set  out  to  determine  such  relationships  by  observation 
and  empiricism. 

With  this  technique,  each  observer  drew  hia  own  conclusions  from  analysi  of 
data  on  either  one  specific  type  and  size  of  fan  or  at  best  a  limited  range  of 
fan  types  and  sizes.  Since  the  relationships  determined  between  r.03.se  and  de¬ 
sign  or  performance  characteristics  were  based  largely  upon  observed  data 
trends,  the  validity  of  those  relationships  was  limited  by  the  ability  of  the 
observer  to  be  aware  of  and  properly  account  for  all  of  the  factors  which 
should  be  observed  for  meaningful  correlations.  The  basic  problem  with  the 
empirical  approach  was  that  it  was  virtually  impossible  to  predict  with 
confidence  the  acoustical  characteristics  of  a  fan  orior  to  building  the 
vehicle  and  obtaining  full-scale  test  data.  Another  problem  was  that  ore  would 
tend  to  assume  that  design  changes  resulting  in  a  noise  change  on  one  f  n  or 
type  of  fan  would  result  in  a  similar  noise  change  on  other  types  of  fa  is. 

Such  assumptions  are,  of  course,  extremely  tenuous,  and  quite  often  erraneous! 

This  program  was  established  in  order  to  develop  meaningful  fan/compressor 
noise  prediction  methods  which  were  based  on  an  understanding  of  fundamental 
noise  generation  and  transmission  mechanisms,  and  verified  by  tests  and 
measurements  on  several  different  scales  and  types  of  vehicles.  The  exper¬ 
iments  were  chosen  not  only  to  establish  functional  relationships  between 
the  basic  mechanisms  and  the  resultant  noise  levels,  but  also  to  evaluate  the 
effects  of  a  series  of  design  parameters  and  concepts  upon  noise  generation 
and  reduction,  and  to  demonstrate  the  accuracy  of  the  prediction  methods. 

This  report  documents  the  results  of  the  program,  starting  with  the  derivation 
of  the  theoretical  and  semi-empirical  relationships  which  describe  the  mechanisms 
of  fan  and  compressor  noise  generation,  transmission,  radiation  and  propagation. 
These  mechanisms  are  described  in  terras  of  aerodynamic  and  geometric  design 
parameters  and  are  quantitatively  related  to  noise  in  terms  of  sound  power  spectra, 
directivity  indices,  and  the  resultant  sound  pressure  spectra. 

The  verification  of  the  validity  of  the  inte^-relationships  between  these 
mechanisms  in  terms  of  pure  tone  and  broadband  spectra  ant  directivities  is 
provided  on  a  number  of  different  fan  and  compressor  vehicles.  Effects  of 
specific  design  parameters,  such  as  blade-vane  numbers,  blade-vane  spacing, 
loading,  tip  speed,  etc,  are  analytical  y  evaluated  and  experimentally  verified. 

Major  subjects  treated  in  the  report  include: 

The,  measured  acoustic  relationship  between  pure  tone  noise  and  blade/vane 
numbers  on  a  controlled  basis. 

Ihe  measured  acoustic  comparison  between  ar  iGV  fan  and  an  IGV-less  fan 
designed  for  essentially  the  same  aerodynamic  performance. 


The  two  and  three  dimensional  cascade  measurements  enabling  relationships 
to  be  derived  between  broadband  fan  noise  and  turbulence  intensities. 

*  The  measured  relationship  between  refraction  phenomena  and  directivity  from 
a  tone  generated  in  a  co-annular  bypass  duct. 

The  measured  acoustic  effects  of  "wake-filling"  relative  to  wakes  generated 
upstream  of  a  rotor. 

*  The  demonstration  of  the  applicability  of  the  basic  mechanisms  over  a  16:1 
size  range. 

*  The  relative  importance  of  blade  loading  and  tip  speed  on  fan  noise. 

"  The  correlation  of  the  effects  of  wave  cancellation  and  reinforcement  pro¬ 
pagation  phenomena  as  a  function  of  microphone  height,  fan  height,  terrain  effects 
and  frequency. 

The  analytical  and  experimental  investigations  in  this  report  result  in 
comprehensive  fan  and  compressor  noise  prediction  systems  which  are  shown  to  be 
reasonably  accurate  and  which  can  provide  a  useful  tool  to  the  engineer  in  designing 
low  noise  fans  and  compressors. 


2 


II.  FAN/COMPRESSOR  NOISE  GENERATING  MECHANISMS 


The  acoustic  signature  of  modern  high  speed  fans  and  compressors  engine  con¬ 
sists  essentially  of  three  different  types  of  noise  components. 

.  Blade  passing  frequency  tones  and  higher  harmo’iics. 

.  Broadband  noise. 

.  Multiple  pure  tones  -  blade  passing  frequency  sub-harmonics. 

This  section  will  describe  the  different  generating  mechanisms  tor  the  three 
noise  components  as  well  as  the  mathematical  models  that  were  developed  to  ex¬ 
plain  and  evaluate  quantitatively  the  fan/compressor  generated  noise  spectra. 


(A)  Blade  Passing  Frequency  Tones  and  Harmonics 


(1)  Mathematical  Model.  The  basis  for  the  analysis  of  pure  tone  genera¬ 
tion  is  a  solution  to  the  linearized  wave  equation  in  an  annular  duct  (see  Figure 
II-A1)  with  uniform  axial  flow  and  solid  body  rotational  fluid  motion.  The 
equation  can  be  expressed  in  r.or  -dimensional  form  as 

2 


9t 


+  M 


9 

9t 


+  M„ 


9 

90, 


(II-A1) 


where 

p  -  sound  pressure 

z  -  axial  distance 

M  -  axial  Mach  number 
z 

M  -  tip  swirl  Mach  number 

v 

9  -  circumferential  angle 

t  -  time 

and  where 


?2 


(II-A2) 


The  solution  to  this  equation  can  be  expressed  as 
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FIGURE  II-A1  DUCT  GEOMETRY  FOR  BLADE  PASSING  FREQUENCY  TONES  AND  HARMONICS  ANALYSIS 


00 


in(e-u)t) 


+  00 

p  *  Z  Z  pc^  noj  e 
n  =  -<o  m15  0 


R  (A  r)  e 
n  tun 


iaz 


[k, 


'z  +  K2  e-1BZ 


where 

„  *”  "2  »9> 

'  1  -  M  2 

z 

b  -V"2  V2  -  a  -  Mz2) 

1  -  M  2 
z 


and  where 


(II-A3) 


(II-A4) 


(II-A5) 


p  -  density 
c  -  speed  of  sound 
n  -  spinning  lobe  number 
<jj  -  mode  tip  Mach  number 
A^  -  eigenvalue 

K1,K2~  constants 

R  (A  r)  -  linear  combination  of  first  and  second  kind  Bessel  functions 
n  run 

The  eigenvalues,  A^j  mav  be  detetuiined  from  the  equation 


Y' 

n 


(A  h)  J'  (*nm)  -  Y’  (A  )  j'  (A  h) 
nm  n  n  nm  n  nm 


(TI-A6) 


where 


h  -  radius  ratio  of  the  annulus 

J*  -  derivative  of  a  first  kind  Bessel  function 

Y*  -  derivative  of  a  second  kind  Bessel  function 
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The  constants  1L  and  K„  are  determined  from  boundary  conditions  at  the  plane 
z  =  0.  These  boundary  conditions  are  based  on  the  annular  duct  geometry  shown 
in  Figure  II-A1. 

It  may  be  assumed  that  at  the  plane  z  *  0  there  exists  a  rotating  line  vortex 
as  shown  in  Figure  II-A2.  It  may  also  be  assumed  that  each  time  the  vertex 
passes  the  angular  position  0  there  will  be  an  Impulsive  change  in  tangential 
velocity.  That  is. 


v  «  f  (r)  6  (0-tot) 


(1I-A7) 


where 


6  (0-tot)  * 


1  when  0  =  tut 
0  when  0  ft  u)t 


and  where  f(r)  is  the  radial  distribution  of  tangential  velocity.  From  equation 
(II-A7)  a  radial  circulation  may  be  computed  as 


T(r)  «*  Rq  cr  f(r)  (II-A8) 

Therefore,  a  boundary  condition  of  the  form 

v/z-0  *  6(9-  wt)  (II-A9) 

z  u  R  rc 
o 

can  be  formulated. 

For  M  such  rotating  line  vortices  equally  spaced,  equation  (II-A9)  can  be  ex¬ 
pressed  as 


v/ 


z=0 


Ihl 

R  c 
o 


M  9  irk 

l  6(0-  tot  +-££-) 

h  =  1 


(II-A10) 


Another  boundary  condition  at  the  plane  z  =  0  is  imposed  by  mass  continuity. 
Since  there  is  no  discontinuous  density  or  area  change  at  z  =  0,  the  axial 
velocity  must  be  continuous.  That  is 

u/z_q  =  continuous  (II-A11) 
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DUCT 


If  a  velocity  potential  is  defined  in  the  annular  duct  as 


♦  ■  e*n  Rn  (*  r)  eiGlz 

Tnm  '  ~n  v  nm 


[h  ^  *  h  ^ 


(II-A12) 


then  the  non-dimensional  velocities  may  be  defined  as 


v  1  3$ 

v  ”  T”  "W 


(II-A13) 


(II-A14) 


Applying  the  boundary  condition  of  equation  (II-A11)  to  equation  (II-A12)  it 
is  found  that 


K1  *  "  K2 


(II-A15) 


Combining  equations  (II-A10), (II-A12),  and  (II-A13)  it  can  be  shown  that 


r (r)M  *  4  i  n  HR  c  A  R  (X  r) 
o  nm  n  nm 


(II-A-16) 


for  each  n,  m  mode.  Multiply  bolt  sides  of  equation  (II-A-16)  by  r  R  (X  r)  and 
integrating  from  h  to  1,  it  is  found  that  n  001 


X  „  -  T  T(r)  r  R  (X  r)dr 

™  U^Rcnl  J.  n  ^ 

o  nm  h 


(II-A17) 


where  n  is  the  norm  of  R  (X  r)  is  the  region  h  to  1. 
nm  n  nm 

Therefore,  the  acoustic  pressure  in  the  duct  may  be  expressed  as 


P  =  £  l  A  P  c2n“>  ein(6"wt)  R  (A  r)  e1®2  e1Bz  (II-A18) 

n  nm  n  nm 

n  ■  -®  tn  =  0 
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The  index  n  is  determined  by  the  well  known  lobe  number  equation 


n  *  n  B  -  kV 


where 

n  -  harmonic  number 
B  -  number  of  blade 
k  -  integer 
V  -  number  of  vanes 

The  index  m  represents  the  radial  mode  structure  within  the  duct. 

The  major  unknown  in  equation  (II-A18)  is  A  which  in  turn  depends  on  T(r) 
in  equation  (II-A17) .  In  the  following  sections  this  unsteady  circulation  is 
defined  for  various  noise  generating  mechanisms. 

(2)  Rotor  Alone  Noise.  The  rotor  alone  noise  is  a  result  of  the 
rotation  of  the  steady  (in  the  rotating  reference  frame)  rotor  pressure  field 
which  is  seen  by  the  stationary  observer  as  a  periodically  varying  (at  the  blade 
passing  frequency)  pressure  pulse. 

An  approximate  calculation  of  the  circulation  about  a  two  dimensional  airfoil 
can  be  obtained  as  follows: 

Consider  the  cascade  of  airfoils  shown  in  Figure  II-A3. 

The  circulation  is  defined  as: 


(II-A19) 
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(3)  Viscous  Wake  Interaction  Noise.  A  viscous  vake  interaction  results 
vhen  the  wakes  from  an  upstream  blade  row  impinge  on  a  downstream  blade  moving 
relative  to  the  upstream  blades.  The  'reaction  of  the  downstream  blade  to  the 
wake  is  a  fluctuation  of  its  loading  (circulation)  which  in  turn  presents  a 
periodically  varying  pressure  source  to  the  stationary  reference  frame.  This 
pressure  fluctuation  is,  therefore,  a  source  of  blade  passing  frequency  noise. 

The  analysis  of  viscous  wake  interaction  generated  pure  tone  noise  is  based 
on  reference  1.  The  analysis  depicts  the  two  blade  rows  as  made  up  of  a  tandem 
cascade  of  thin  airfoils  in  an  inviscid  two  dimensional  shear  flow  with  the  wake 
shape  empirically  determined. 

Figures  II-A4  and  II-A5  contain  the  basic  geometric  configuration  used.  The 
wake  shape  can  be  represented  by 


Y  =  .962  C*/2  c  (-  -  .7) 1/2  (II-A22) 

d  c 


where 

Y  -  wake  semi-width 
c  -  blade  semi-chord 
-  drag  coefficient 

x  -  distance  downstream  from  the  upstream 
blade  row  mid-chord 

Using  the  geometry  of  Figure  II-A5  and  taking  the  relative  motion  of  the  two 
blade  rows  into  account  the  unsteady  upwash  on  the  rotor  can  be  expressed  as 


is  ** 

V  =  tr  Vr  -  1  G-  exP  ivrn  (t  -  v  )  (II-A23) 

n  =  1  n  L  r  . 


where 


2nU 


s 


(II-A24) 
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FIGURE  II-A4  BLADE  WAKE  PROFILE 


g! 


v 

14.622  ^ 


r 


C.  sin  0o 
d  s 

Or  +  -3)  cos  a, 

s 


(II-A25) 


c 

s 


sec  og 


(II-A26) 


(II-A27) 


and  where  o.  is  the  stator  solidity  (all  other  variables  are  defined  in  Figure 
II-A5) 

It  might  be  interesting  to  note  that  equation  II-A27  implies  that  the  wake  widtl 
and  velocity  uc  are  measured  at  the  quarter  chord  point  on  the  rotor.  Variations 
in  this  point,  however  do  not  affect  the  results  significantly. 

Using  the  results  of  thin  airfoil  theory,  (reference  2)  the  unsteady  circulatioi 
can  be  obtained  from  the  upwash  velocity  as 


r| 


c  V 
r  r 


j(3^|  1 J (n  w)  |  |s(n  «>)  i 


(II-A28) 


|j(nu))  | 


(n  ui) 


+ 


Ji 


1/2 


(II-A29) 
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ll)  3 


c  iU 

r  s 


(II-A30) 


c 

3 


and  where  | S  (n  w) ;  can  be  approximated  by 


S 


(nw)  | 


a  +  nm 

(tra  +  1)  nui  +  2nn2(u2 


1/2 


(II-A31) 


where 


a  =  .1811 


The  circulation  calculated  by  equation  (II-A28)  is  the  required  unsteady  cir¬ 
culation  due  to  the  interaction  of  the  wake  of  an  upstream  stator  with  a  downstream 
rotor  at  a  given  radius.  Similar  equations  can  be  derived  for  the  interaction  of 
a  rotor  and  downstream  stator  by  interchanging  the  subscripts  r  and  s  in  equations 
IIA-23  through  IIA-30. 

In  order  to  get  an  approximate  radial  distribution  of  unsteady  circulation,  the 
calculation  can  be  made  at  three  radial  stations:  hub,  pitch  and  tip,  and  a  curve 
fit  through  the  results.  Then  the  radial  distribution  can  be  expressed  as 


r(r)  —  a  +  a. r  +  a0r2 

o  1  2 


(II-A32) 


In  sumnu.;~y,  the  unsteady  circulation  has  been  cast  in  such  a  form  that  it  can 
be  directly  related  to  the  parameters  which  describe  the  geometric  and  aerodynamic 
configuration  of  the  blade  rows. 

Equation  (II-A32)  can  be  substituted  in  equation  (II-A17)  to  obtain  the  ampli¬ 
tude  coefficient  A  . 

tun 
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(4)  Potential  Interaction  Noise.  A  potential  Interaction  between  blade 
rows  produces  a  mutual  disturbance  and  thus  two  noise  sources  are  generated. 

That  is,  the  pressure  field  of  the  rotor  disturbs  the  stator  and  the  pressure 
field  of  the  stator  disturbs  the  rotor.  To  a  good  approximation  these  two  noise 
sources  can  be  thought  of  as  independent.  These  disturbances  also  lead  to 
periodic  blade  loading  (circulation)  variations  which  produce  pure  tone  acoustic 
waves . 

Two  analyses  were  developed  for  potential  interaction.  One  two  dimensional 
and  the  other  a:.isymmetric. 

The  two  dimensional  analysis  was  adapted  from  reference  2.  The  intent  of 
this  work  was  to  determine  the  upwash  (or  downwash)  on  a  blade  due  to  the  poten¬ 
tial  field  of  a  neighboring  blade  row  moving  relative  to  the  blade.  The  resulting 
upwash  velocity  can  then  be  used  in  a  calculation  of  the  unsteady  circulation. 

The  analysis  is  for  a  two  dimensional  cascade  of  thin  airfoils  in  a  non-vlscous 
incompressible  flow.  An  approximate  three  dimensional  solution  can  be  obtained 
by  making  a  two  dimensional  calculation  at  several  radial  locations  and  fitting 
a  curve  through  the  results. 

Consider  the  two  cascades  shown  in  Figure  II-A6.  The  upwash  velocity  on  an 
upstream  blade  due  to  a  point  vortex  of  strength  r'  on  each  downstream  blade  can 
be  expressed  as 


v 


N 

l 

J  »  1 


(II-A33) 


The  point  vortex  can,  l-wever,  be  replaced  by  a  line  vortex  distributed  along 
the  downstream  blade  row  chord  where 


r 


r 


Y  (x  )  dx 
r  r  r 


(II-A3  ) 


and  y  is  the  distributed  vortex  strength.  This  strenjhc  in  turn  may  be  expressed 
in  terms  of  a  flat  plate  or  elliptic  distribution  along  the  chord.  That  is 


or 


2V  Ar 
r  o 


1  -  cos  6 

A 

sin  6 


y  *»  4V  A,  sin  6 
'r  r  1 


(flat  plate) 


II-A35) 


(elliptical) 


(II-A36) 
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where  the  transformation 


x  «  c  cos  0 
r  r 


(II-A37) 


X  X 

has  been  made  and  where  A  and  A.  are  determined  by  the  blade  loading. 

01 

Substituting  equations  (II-A34-37)  into  equation  (II-A33),  introducing  the 
geometry  of  Figures  II-A6  and  integration  over  the  downstream  blade  row  chord, 
the  upwash  velocity  can  be  expressed  as 


r  e-la8 


1  +  2  l  H-  exp 
n  * 


r  ' 


b  +  ih  -X  e  *as 
.  s  s  > 


(II-A38) 


where 


H-  =  J  (z)  +  i  J  (z  ) 
n  or  i  r 


(II-A39) 


for  a  flat  plate  vorticity  distribution  and 


„r  2iei0lr  ,  , 

Hn  "  Rc"  "  J1  <zr) 


(II-A40) 


for  an  elliptic  distribution.  The  argument  is  expressed  as 


z  -  —  i  [tt /2  +  a  ] 

r  =  too  e  r 

r 


(II-A41) 
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Thus  far  the  airfoils  have  been  assumed  to  be  stationary  relative  to  one 
another.  Introducing  moving  coordinates  the  upwash  can  be  expressed  as 


s  rre  8  .  r,  :  „s  2,  5  -is 

- - 53-  +  E  Gj  exp  -3-  e  s  X; 

r  s  n  \  r 


(II-A42) 


where 


n 


ira  d  . 

— •— — -  e  as  H-  exp  -  tfna 
d  n 

r 


rj^c”  +  ^  tanar  “ 


i£ _ \ 

AA.  rv  ^ 


V  cos  a 
s  s 


-Ml 


(II-A43) 


It  can  be  shown  (references  3  and  4)  that  the  magnitude  of  the  unsteady 
circulation  for  sinusoidal  motion  is  given  by 


r  -  rr  I  |g^|  |j(Xgn)|  | S(nws) 


(II-A44) 


where 


d 


e  [tt/2  -  ag] 


(II-A45) 


2ttU  c 

5 

o  <■  — ;r-, — ; 
s  V  d 
s  r 


(II-A46) 
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|j(Xgn)|  «  |jo(XgH)  -  i  Jx(Xsn) | 


(II-A47) 


|S(nw  )|  ■  - 
3 


a  +  nw 

t 


s 


_T 


,1/2 


-  j. 


(II-A48) 


and  where 


a  «  .1811 


(II-A49) 


Therefore,  the  unsteady  circulation  can  be  calculated  from  equation  (II-A44) . 
For  an  axisymnetric  cascade,  the  circulation  can  be  obtained  approximately  by 
applying  equation  (IX-A44)  at  several  radial  locations.  The  resulting  circulation 
distribution  can  then  be  applied  to  equation  (1I-A17)  to  determine  A^. 

Table  II-A1  Indicates  the  changes  to  be  made  in  equations  (I1-A39-46)  to  ob¬ 
tain  the  remaining  three  types  of  potential  interactions. 

For  the^axisymmetric  analysis  consider  an  element  dR  of  a  vortex  line  of 
strength  r  located  at  a  point  on  the  rotor.  Due  torthis  vortex  element, 
the  induced  velocity  at  a  point  Rs  on  the  stator  is  given  by  (Figure  H-A7) 


r  •  x  p 
r  r  r 
4irp3 


(II-A50) 


or  in  scalar  form 


d<U 


r'(x.  R.)  sin6 

- Lj- - 


fefr2 


dR 


(II-A51) 


The  induced  velocity  component  vertical  to  the  chord  of  the  stator  blade 
Figure  (II-A8)  is  then 


du 

s 


rr  (xr>Rr)sin(S 


4w 


2 


sinn  dR 

r 


(II-A52) 
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SUBSCRIPTS 

r  or  s 

FLOW  ANGLES 

a 

WHEEL  SPEED 
U 

IGV  -  Rotor 
Potential 

(as  shown) 

(as  shown) 

(as  shown) 

Rotor  -  IGV 
Potential 

r  -*•  s 

a  ■+  -a 

U  ->  -U 

OGV  -  Rotor 
Potential 

(as  shown) 

a  -*•  -a 

u  -u 

Rotor  -  OGV 
Potential 

r  -*•  s 

(as  shown) 

(as  shown) 

TABLE  II-A1 


Subscript  Changes  for  Multiple  Potential 
Interactions 
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STATOR 


IGURE  II-A7  AXIAL  VECTOR  DIMENSIONS  FOR  AN  AXISYMMETRIC  POTENTIAL 
INTERACTION 


where 


sinS  -  %  J[Z2  +  (R  -  R  )2]1/2 
s  s  r  s 


(II-A53) 


Z 

s 


+  (Z  J 
s  2 


(II-A54) 


(*),  *  (C  +  x  )  cos  a  +  (C  _  x  )cosas  +  b' 
s  x  r  r  r  s  s 


(II-A55) 


(Z  )0  *=>  (C  +  x  )  sinar  -  (C  -  x  )  sinas  +  ndr  -  Ut  (II-A56) 
s  L  rr  *•  ss  *• 


Z  is  the  distance  between  a  point  on  the  stator  and  a  point  on  the  rotor.  U  is 
tBe  speed  of  the  rotor,  positive  down  (Figure  II-A8). 

The  following  relations  will  be  used  in  the  later  development: 


V  sin  (a  +a  )  =  U  cosar 

s  r  s  r 

V  sin  (a  +  a  )  *  U  cosa„ 

r  r  s  s 


(II-A57) 


The  angle  n  between  the  chora  of  the  stator  blade  and  the  velocity  qs  can  be 
expressed  in  terms  of  ag,  an-  '^s^2» 


n  = 


[a  +  tan-1  ^^2  1 

L s  tv,  J 
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BY  THE  STEADY 


now 


sin i)  *  { (C  +  x  )  cos  (a  +  a  )  +  b’  cosas  +  C  -  x  -nd  sinas 

r  r  r  8  s  s  r 


+  Ut  sin  a,  }/£ 


(II-A58) 


The  vertical  velocLty  at  the  stator  induced  by  the  steady  non-uniform  cir¬ 
culation  distributed  along  the  radial  direction  of  the  rotor  blade,  is  obtained 
by  integration  equation  (II-A52)  from  R r  =  to  R£. 


u 

s 


JL 

4ir 


^xr»Rr^  sin6 


(-la.  \2 
^sin6 


sinn  dRr 


(II-A59) 


Substitution  leads  to 


u  s 

s 


(xr’Rl^ 


(£_). coso  -  (£  )2  sin  a 
si  b  0  « 

uF 


+  (r  -  R  )zrj/'J 

r  o 


dRr  (II-A60) 


If,  now,  T  is  replaced  by  the  distribution  of  vortices  y (x  ,Rf)  dx 
along  a  line  representing  a  typical  rotor  blade,  the  induced  vertical  velocity 
at  the  stator  due  to  a  single  rotor  blade  is  obtained  by  integrating  over  x^_ 
from  -Cf  to  Cf. 


u 

s 


Y  (x  ,R  ) 
r  t  r 


(£g)l  cosas  -  (£g)2  sinag 

[l  '  +  (R  -  R  dRrdxr 


(II-A61) 


By  assuming  for  the  present  computation 


y  =  (a,  +  a0R  +  a,R  O 
r  1  2  r  3  r 


2 \  1  -  cos  0 


* —  r 

sin 


(II-A62) 
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where  the  relationship  between  *r  and  0r  Is  given  by 


xr  *  cos  ®r 


(II-A63) 


Equation  (II-A61)  can  be  written  as 


ir  R, 


(1  )i  C08as  -  (£  >2  as 


1  "  '  o'*  - o  '  cr  *  M 

-  f  f  (^i  +  +  aJt  2)  /  • 

*  Jo  h  Hs2  +  (Rr  -  Rs>2  13/2 


(1  -  cos0r)  CrdRrd0r 


(II-A64) 


Summing  up  the  contributions  of  all  rotor  blades ,  the  induced  vertical  velocity 
at  the  stator  blade  due  tc  the  whole  row  of  rotor  blades,  is  given  by 


n?  t;  R 

u  =  1  l  f  f  2  (a  +  a2R. 


r  +  a3  Rr2)  * 


r  1 


(1  )i  cosas  -  (£s^2  sin«8 

— § - 7 -  (1  -  cos0r)  C  dR  de 

[I  2  +  (P  -  R  )2)3/2 
s  s 


(II-A65) 


where 


B  -  1 
n,  =  -  _r _ 

1  2 


B  -  1 

n  «  — — ^ —  if  B  is  an  odd  number 


B  -  2 

-  r 


n2  2 


if  B  is  an  even  number 
r 
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Observing  that  (1  ) i ,  (t  >2,  1  are  quantities  evaluated  at  Rr  =  Rg,  the 
Integration  can  Se  carried  out  analytically.  r  S 


n2  3 


u  (0  ,R  ,t)  =  -7-  1  /  H  (0  ,0  ,  R  , 

s  s  s  __  4  r  r  s  s 


t)  lUs)i  cos  as  -(Jts)2sinasl 


n=n^  r=0 


(1  -  cos  0_)  C  c 
r  r 


(II-A66) 


where 


x  C  cos  0£ 
s  s 


.  .  R2-  R  Rl  -  R  "I  a2  Tr2R  -  R  2  -£  2 

H  (6  ,0  ,R  ,t)  -  fV  - £  _  - £  +  —7  —2 - § - a 

r  r*  s’  s’  l ^  L  p2  pl  J  z3  |_  P2 


RtR  -  R  2-  l  2 
1  s  3 _ s 


1  U-F* 

_  3  j  P2  _ 


.  2-  l  2 

Vr^  (r2  -  V 


-2R 

SJ 


r  R  2-  l  2  1  R2  _  R  +  p2  ) 

-  TT  -V-  (R>  -  V  -  2RS  +  1orIr,  -R%pJ  (II-A67) 

L  s  J  s  ' 


^s+(Rl  -Rs)2 


(II-A68) 


/  2  2 

/*  +  (R2  -  R  ) 

s  s 


(II-A69) 
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In  order  to  evaluate  the  circulation  of  the  stator  blade  due  to  the  effect 
of  the  rotor,  equation  (I1-A66)  nay  be  written  in  the  form: 


■  vs 


„  unv  1 1 
C  e  si 

m  l 


(II-A70) 


where  £  '  indicates  the  summation  over  m  omitting  m  =  0  term, 
m 


U 

s 


V 

s 


-imv-t 
e  8 


dt  ,  p  =  ir/v^. 


vs  =  2it0/dr 


The  time-independent  term  in  u8  of  equation  (II-A70)  represents  the  steady 
induction  effect  of  the  rotor  on  a  stator  blade.  It  will  be  omitted  in  the 
later  development  for  only  the  unsteady  effects  are  of  Interest  here. 


The  coefficients  C  which  are  functions  of  Ogycan  be  represented  by  a  cosine 
series:  m 


C 

m 


=  A  +2 


mo 


cos  q6g 


(II-A71) 


where 


e-imvst 


cos  qe 


dtde 

s 


Substitution  of  equation  (II-A71)  into  equation  (II-A70),  (omitting  the  term 

C  ) ,  leads  to 
o 


u 

s 


V 

s 


00 


00 


V 


m= 


£ 

q=l 


cos 


imv_t 
e  s 


(II-A72) 
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Sears  (reference  4)  has  shown  that  when  the  relative  vertical  velocity  u(x)  of 
the  airfoil  is  given  in  the  form 


k 


+  2  l  Aq  cos  qe  VVe 

q=l  ‘ 


ivt 


the  lift  is  given  by 


(II-A73) 


L  *  2irpCV2  eivt  |(Aq  +  A^  C  (u»)  +  (Aq  ”  A2>  ^  }  (H-A74) 

where 

x  =  C  cos  0 

V  =  the  incoming  air  velocity 
p  =  the  density  of  air 
C  *  the  serai-chord  of  the  airfoil 
u)  *  the  reduced  frequency  *  v  c/V 

v  »  the  frequency  of  oscillation  in  radians  per  second 
Kl(la)) 

C(a>)  =  K  (iui)  +  K.  (iuj) 
o  i 


In  order  to  obtain  the  quasi-steady  circulation,  the  results  of  Von  Karman 
and  Sears  (reference  7)  will  be  utilized.  It  states  that,  if  the  motion  of  the 
airfoil  is  such  that  its  quasi-steady  circulation  is  given  by 


rq  (t)  =  F  eivt  (II-A75) 

then 

Lq  +  L(2)  =  Lq  C (to)  (II-A76) 
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where  Lq  is  the  quasi-steady  lift  and  is  the  lift  due  to  wake  effects. 
Comparison  of  equation  (II-A74)  with  equation  (II-A76)  reveals 


Lq  -  2irpc  V2  e  lvt(Ao  +  kj  (II-A'7) 

Since 

Lq  -  prq  V  (1I-A78) 

the  quasi-steady  circulation  is  given  by 

rq  ®  2nC  V  e  lvt  (Ao  +  A1)  (II-A79) 


As  mentioned  in  reference  6f  using  the  results  of  Von  Karman  and  Sears 
(reference  7),  it  is  easy  to  show  that  the  total  circulation  is  given  by 


r(t)  -  rq  (t)  e"lu  S(ut) 


(I1-A80) 


Comparison  of  equation  (II-A72)  with  equation  (II-A73)  and  using  equation 
(II-A80)  gives 


r8  *  2*  Csvs  1  {eln(Vst-“s>  s 


m«l 


+  e 


-im(vst  -  ug) 


S  (~mu  )  (A 

s  -mo 


(II-A81) 


The  velocity  at  the  rotor  induced  by  the  steady  circulation  of  the  stator  can 
be  obtained  in  a  similar  manner  as  that  <ised  to  determine  the  induced  velocity 
at  the  stator,  Figure  1I-A9. 
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The  result  is 


r  =  2ttC  V  L-'{eim(vrt~wr)  S  (sm  )  (B  +  B  )  +  e"im(vrt'wr)  . 
r  r  r  v  r  mo 


ra=l 


3  r>  «-o  '  B-ml>  } 


(II-A82) 


where  B  *s  are  the  coefficients  of 
mq 


u  =  V 
r  r 


1  1 
m=  — » 


B  +  2  Z 
mo 

q*i 


B  cos  q0 
mq 


4 


elmV 


(II-A83) 


e 


-imv^t 


cos  qerdtde„ 


(II-A84) 


vr=  2  nU/ ds 

oir  =  vrCr/Vr 


Application  of  equations  (II-A81)or  (II-A84)  at  several  radial  stations  along 

the  rotor's  span  will  provide  a  distribution  of  unsteady  circulation  which 

can  be  used  to  evaluate  A  , 

ran 
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The  application  of  the  analytical  results  indicate  that  in  most  cases  the 
two  dimensional  and  the  axisymmetrlc  analyses  yield  comparable  data.  It  is 
only  for  very  closely  spaced  rotors  and  stators  (less  than  0,1  chord  spacing) 
that  the  axisymmetrlc  analysis  is  required. 

(5)  Wake  Vortex  Interaction  Noise.  A  secondary  effect  of  the  blade 
row  interactions  is  a  periodic  wake  vortex  interaction.  A  schematic  of  this 
phenomenon  is  shown  in  Figure  (II-A10).  The  shed  vortex  is  usua.^v  termed  a 
"starting  vortex"  and  is  a  necessary  result  of  the  conservation  of  vorticity 
in  a  closed  conservative  system.  Thus,  when  two  blade  rows  interact  two 
starting  vortices  are  shed  into  their  wakes.  The  vortices  in  turn  impinge 
on  the  downstream  blade  row  and  set-up  still  another  circulation  fluctuation. 

It  can  be  shown  (reference  2)  that  the  wake  vortex  strength  due  to  the  poten¬ 
tial  interaction  of  a  rotor  pressure  field  with  an  upstream  stator  is  given  by 


s 


e. 


co  — 

t)  =  -  rr/cg  _Z  K(nws)J(nAs)  e  nys 
n*l 


"  Xs/Vs) 


(II-A85) 


where  G-  is  given  by  equation  (II-A43)  and  x°  is  a  distance  measured  along  the 
stator  chord. 

This  strength  may  then  be  translated  i..to  an  effective  induced  velocity  at 
the  leading  edge  of  the  rotor.  The  induced  velocity  can  be  expressed  as 
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CO 
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n  = 


CO 
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1  k=  -  “ 
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nk 


o 

x 

exp  [  iyrk  (t  -  ■—)] 
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(II-A86) 


where 


w-, 

nk 


r  G-  K  (no)  )  J  (nAs)  exp 
r  ° 


[-innaj 


(bU/cgVg  cos. 


-)  +  it 
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doCosa,- 


(Vsdr/Vr'is)  (k/n) 


1  +  cot?  $[1  -  (V  dr/V^.dc)  (k/n)  sec  3] 


TiT 


r  s' 


(II-A87) 


It  may  be  noted  that  the  sum  over  k  excludes  the  tern  k  =  0. 
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A  POTENTIAL  INTERACTION  TAKES  PLACE 


STATOR  ROTOR 


STATOR 


A  VORTEX  IS  SHED  FROM  THE  UPSTREAM  STATOR 
STATOR  ROTOR  STATOR 


VORTEX  CHANGES  THE  DOWNWASH  VELOCITY  ON 
DOWNS  TRiAM  ROTOR 


FIGURE  II-A10  SCHEMATIC  OF  A  WAKE  VORTEX  INTERACTION 


'<4 


Using  the  results  of  thin  airfoil  theory,  the  unsteady  induced  velocity  may 
be  expressed  in  terms  of  an  unsteady  circulation  magnitude. 


co  oo 

|r|  =  j:  l  2ncr|w-k|  |j(iok)|  |s(u)_k)|  (II-A88) 

n=l  k=  -»  u 


Equation  (II-A88)  can  then  be  applied  at  the  hub,  midspan,  and  tip  of  the 

rotor  to  obtain  a  radial  unsteady  circulation  distribution  which  in  turn  may 

be  used  to  calculate  A  in  equation  (II-A17) 

nm 

A  similar  result  can  be  derived  for  an  upstream  rotor-downstream  stator  com¬ 
bination.  It  is  only  necessary  to  reverse  the  r  and  s  subscripts  in  all  the 
equations . 

(6)  Summary .  Table  II-A2  summarizes  the  information  required  to  make 
a  calculation  of  blade  passing  frequency  and  harmonic  no-^se. 

Figure  (II-A11)  shows  a  typical  calculated  result  of  spacing  versus  noise  for 
the  four  types  of  interaction  noise  possible  for  a  rotor  and  OGV.  In  this  case 
viscous  wake  interaction  is  dominant  for  about  .15  chord  spacing  or  greater.  The 
viscous  wake  interaction  also  shows  the  slowest  decay  rate  with  spacing. 

It  must  be  pointed  out  that  different  aerodynamic  designs  will  have  different 
spacing  decay  rates.  The  parameters  listed  in  Table  II-A2  and  the  equation  for 
the  unsteady  circulation  clearly  show  the  importance  of  the  details  of  the  aero¬ 
dynamic  design  on  the  tone  noise  generation.  The  mathematical  models  described 
above  were  very  successful  in  providing  a  good  quantitative  prediction  of  blade 
passing  frequency  tones  and  their  higher  harmonics.  Detailed  comparisons  between 
prediction  and  experimental  data  obtained  on  a  number  of  vehicles  are  presented 
in  Section  VIII. 
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*blade  row 
solidities 

tip  radius 

hub  radius 

number  of  blades 

number  of  vanes 


rotor  speed 

speed  of  sound 
density 

*relative  and  absolute 
air  angles  entering  and 
leaving  each  blade  row 

*sectlon  drag  coefficient 


number  of  spinning 
lobes 


harmonic  number 
(1  =  fundamental) 

axial  distance  from 
t.ha  plane  of  generation 

*axdal  spacing  between 
blade  rows 

*  required  at  hub,  midspan  and  tip  radii 

TABLE  II-A2 

Data  Required  for  a  Blade  Passing  Frequency 
PWL  Prediction 
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VISCOUS  WAKE  INTERACTION 


xarao. 
ini  ova 


NOMENCLATURE  SECTION  II-A 


A  coefficients  of  relative  vertical  velocity  distribution  - 

equation  (22) 

A^  coefficients  of  the  chordwise  vorticity  distribution 

A  n.  m  mode  amplitude  constant 

nm 

b  blade  row  axial  spacing,  mid-chord  to  mid-chord 

b'  blade  row  axial  spacing,  trailing  edge  to  leading  edge 

c  speed  of  sound  or  airfoil  semi-chord 

C^  drag  coefficient 

Cu  steady  tangential  velocity 

C(w)  Theodorsen  function 

d  blade  to  blade  spacing 

e  natural  log  base 

f(r)  radial  distribution  of  tangential  perturbation  velocity 

G  coefficient  of  unsteady  wash 

n 

h  radius  ratio  of  the  annulus 

hg  see  Figure  II-A6 

H-  see  equations  (II-A39)  and  (II-A40) 

n 

i  imaginary  number  indicator 

J  see  equation  (II-A47) 

first  kind  Bessel  function  of  first  order 

J  first  kind  Bessel  function  of  zeroth  order 

o 

constants 

£  distance  between  a  point  on  the  stator  to  a  point  on  the  rotor 

at  the  same  radial  distance 
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L  lift 


m 

M 

M 

2 

M 

o 

n 

n 


R 

R, 


n  (Anm 


R 

r 

S 

t 


r) 


u 

u 

c 

U 


V 

V 
X 

y 

Y 

Y 

n 

z 

r 

Z 


radial  pressure  mode  index 

number  of  rotating  line  sources 

axial  Mach  number 

tip  swirl  Mach  number 

spinning  lobe  number 

harmonic  number 

pressure 

radial  distance 

radial  dimension 

hub  radius 

tip  radius 

linear  combination  of  first  and  second  kind  Bessel  function 

see  equation  (1I-A31) 

time 

axial  perturbation  velocity 
wake  velocity  decrement 
linear  rotor  speed 

tangential  perturbation  velocity  or  upwash  velocity 
free  stream  steady  velocity 

distance  measured  parallel  to  the  blade  chord 
distance  measured  perpendicular  to  the  wake 
tvake  semi-width 

second  kind  Bessel  function  of  order  n 
see  equation  (II-A41) 
axial  distance 
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a  - 
0  - 
0  = 
r  - 
r*  - 

Yr  - 
6  - 
£  - 

4  - 
n  - 

n 

nm 

e  _ 


3 


nm 

A 

s 

u 

C 

P 


P1P2 


a 


r 

n  _ 

0)  _ 


angle 

see  equation  (II-A5) 

a  +  a 

s  r 

circulation 

point  vortex  circulation 

vortex  distribution  along  the  rotor  chord 

Krondf.cker  delta  function 

see  equation  (II-A85) 

complex  distance  from  rotor  to  stator 

distance  measured  tangentially  along  the  blade  rows 

norm  of  the  cylinder  function 

circumferential  angle 

transformed  coordinate  along  the  airfoil  chords 
eigenvalue 

see  equation  (II-A45) 

circular  frequency  of  motion  in  time  =  WU/d 

dxial  distance  from  a  point  on  the  rotor  to  a  point  on  the  stator 
density 

see  equationg  (II-A68)  and  (II-A69) 
solidity 

velocity  potential 
circular  arc  camber 
rotational  speed  of  rotor  in  rpm 
mode  tip  Mach  number  or  reduced  frequency 
sta  ir  solidity 
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SUBSCRIPTS 


r  rotor 

s  stator 


SUPERSCRIPTS 

r  rotor 

s  stator 

q  quasi-steady 

(2)  wake  effects 
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(B)  Analytical  Prediction  of  Broadband  Noise.  Modern  aircraft  engines  are 
characterized  by  numerous  noise  sources.  One  category  of  noise  which  is  of  con¬ 
siderable  importance,  especially  for  noise  evaluation  of  the  high  by-pass 
tur'oofan  engine,  is  broadband  noise.  The  name  "broadband"  noise  is  derived 
because  of  the  broad  frequency  range  associated  with  the  sound  pressure  level 
spectra  of  aircraft  fans  and  axial  compressors. 

It  is  known  that  there  are  two  significant  sources  of  broadband  noise  in 
fans.  One  arises  in  any  type  of  flow  from  the  vortex  shedding  at  the  airfoil 
trailing  edge  which  induces  local  surface  pressure  fluctuations  on  the  blade. 

The  other  is  the  broadband  noise  generated  when  the  blade  moves  in  a  turbulent 
airflow.  Under  inlet  turbulent  flow  conditions,  the  turbulent  velocity 
fluctuations  ahead  of  the  blade  cause  random  changes  of  incidence  of  the  re¬ 
lative  flow  and  hence  a  randomly  fluctuating  lift  on  the  blade. 

The  two  broadband  noise  mechanisms  have  the  basic  common  feature  of  being 
generated  by  a  random  noise  source  distribution  on  and  in  the  vicinity  of 
a  rotating  blade.  The  basic  acoustic  source  model  associated  with  the  two 
types  of  mechanisms  discussed  above  is  the  dipole  noise  source.  In  addition 
to  a  dipole  random  distribution  which  results  in  a  broadband  noise  spectra, 
there  is  also  a  random  quadrupole  contribution  resulting  from  the  turbulence 
presence  in  a  wake.  The  following  discusses  a  number  of  theories  based  on  the 
noise  mechanisms  discussed  above.  Two  main  catagories  will  be  considered. 
Broadband  formulations  for  predicting  overall  sound  powers  and  formulations  for 
predicting  spectral  sound  level  contributions  of  broadband  noise. 

(1)  General  Formulation.  In  a  completely  general  way  the  wave  equation 
for  broadband  noise  can  be  obtained  from  the  following  conservation  equations: 

Continuity  Equation 


_i£ 

at 


+  Div  ou  =  0 


(II-B1) 


Momentum  Equation 


9pu 

at 


+  Div  (pu)  u  =  Div  x 


(II-B2) 


where  r  <  -p  I  +  t ' 

t'  :  Viscous  Stress  Tensor  [udef  u  +  Adiv  u] 


4  3 


y  *»  First  viscosity  coefficient 
X  “  Second  viscosity  coefficient  (X  =  2/3y  ) 

I  ■  Unity  tensor 

def  u  =  Deformation  tensor  (grad  u  +  (grad  u)^> 

By  using  the  time  derivative  of  equation  (II-Bl)  and  the  divergence  of 
equation  (II-B2)  and  elimination  of  the  momentum  flux  term  pu  one  obtains 


>  > 

=  Div  Div  [  r*  -  (pu)  u] 


(1I-B3) 


If  p  and  p  are  considered  as  constant  reference  values  such  that  re¬ 
ference  speed  of  sound  Cq  _  Yp  /p  *  then  equation  (II-B3)  may  be  written  as 


V  (P  "  p0)  "  C 


=  Div  Div  [t?  -  (pu)  u] 


-  TTT  (p  -  C  2p) 

C  1  3^  v  o  M 

o  L 

This  relationship  represents  the  inhomogeneous  wave  equation  for  pressure 
fluctuations.  Lighthill  (reference  1)  has  shown  that  except  perhaps  in  ex¬ 
ceptional  circumstances,  the  source  term  -  1  92p  (p  -  C2p)  is  negligible. 

Thus,  to  this  approximation  one  has  CG 2  Sp2 


V2p  -  ~2  |-]r  =  Div  Div  (T*  -  (Pu)  u]  (II-B4) 

o  3t 


The  solution  of  equation  (II-B4)  may  be  obtained  in  a  closed  form  by  the 
method  of  "Green's  Function"  (reference  2).  Let  G(r,  t/r  , t  )  be  a  solution 
to  J  ° 


v2g  -  T?  =  "  (>r  -  >ro)  6  (t  -  t0) 

0  t 
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where  6(r  -  r  )  is  the  three-dimensional  apaclal  Dirac-delta  function  and 

6  (t  -  t  )  is  the  time  dependent  form  of  the  Dirac-delta  function, 
o 

Then,  if  p  (r,t)  is  a  solutiojg  gf  the  following  non-hotnogeneous  wave 
equation  with  a  source  function  p(r,t) 


V^p 


(r,t) - —  Y7  p  =  -4trp  (r,t) 


for  t  >0  within  and  on  the  surface  S  as  shown  in  Figure  (II-B1)  one  has 
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(II-B5) 


where  the  subscript,  o,  on  the  vector  "del"  operator  specifies  differentiation 
with  respect  to  the  source  position  coordinates. 

'V  > 

Since  the  source  term,  P(r  ,t)>  and  the  surface  effects  are  of  primary  interest, 
the  volume  integral  expressing  initial  effects  will  be  considered  negligible. 
Equation  (II-B5)  becomes  then 
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(II-B6) 


> 

where  n.  is  the  unit  normal  vector  pointing  outward  from  the  surface  S.  In  thee- 
dimensions,  the  Green's  function  is 


G(r,t/r,tJ  ="  S[R/C  -  (t-t^)];  R2  - (x-xq)2  +  (y-yo)2  +  (z-zo)2 


o  o 


Therefore  for  the  source  term  as  given  in  equation  (II-B4)  one  has 
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(II-B7) 


Integrating  over  time*  yields 


*Brackets  about  terms  signify  the  term  is  evaluated  at  a  retarded  time  (t-R/C) 
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For  far  field  terms: 


Div  Div  [t*  -  (pu)  u)]  dV 
00  o 


Noise  resulting  from  fluctuating  shearing  stresses 
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Noise  resulting  from  the  effect  of  solid  boundaries  upon  the  flow 


(II-B9) 
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Thus  equation  (II-B9)  represents  the  acoustic  pressure  in  the  far  field 
measured  by  an  observer  at  r,t  due  to  fluctuating  shearing  stresses  and  due 
to  the  reactive  effect  of  solid  boundaries  upon  the  fluid  at  r^.t^. 

The  above  formulation  is  derived  from  Curie,  who  performed  the  first  theo¬ 
retical  work  dealing  with  sound  radiation  from  rigid  boundaries.  Based  on 
dimensional  arguments  Curie  had  shown,  that  for  surfaces  small  compared  with 
a  wave  length,  the  total  dipole  acoustic  power  radiating  from  a  solid  surface 
in  a  flow  was: 


l^TTpC2 


i£ 

at 


(II-B10) 


Based  on  equation  (II-B10),  quantitative  estimates  of  the  output  from  a  fan 
blade  can  be  made  by  considering  the  mechanisms  involved  and  by  properly  evalua¬ 
ting  the  integral. 

(2)  Calculation  of  Broadband  Noise  Overall  Acoustic  Power.  In  order 
to  evaluate  equation  (II-B10) ,  Sharland  (reference  3)  performed  dimension  ar¬ 
guments  which  lead  to  expressions  for  the  total  acoustic  power  for  noise  genera¬ 
ted  from  blades  operating  in  laminar  and  turbulent  inlet  flow.  The  following 
paragraphs  discuss  Sharland* s  approach  and  the  application  of  the  analysis  to 
engine  noise  calculations. 

(a)  Noise  from  a  Blade  Operating  in  Laminar  Flow.  Implicit  in  the 
derivation  of  equation  (II-B10)  is  the  assumption  that  it  is  normal  pressure 
fluctuation  on  the  surface  which  are  dominant.  The  fluctuations  due  to  tan¬ 
gential  stresses  are  assumed  to  be  negligible  in  comparison  to  the  normal  pressure 
fluctuations.  Thus  the  fluctuating  pressure  difference  can  be  thought  of  as  a 
local  lift  fluctuation  per  unit  area.  Writing  the  pressure  distribution  in  terms 
of  a.  local  lift  acting  on  a  two-dimensional  section  of  the  blade  one  has 
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where  U  is  the  local  mean  velocity  parallel  to  the  surface  of  consideration. 

Then 
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Upon  substitution  into  equation  (1I-B10)  and  integrating  in  the  streamwise 
direction,  the  expression  for  the  total  acoustic  power  becomes 
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(II-B11) 


where  C  is  the  blade  chord  at  spanwise  coordinate  y,  is  the  fluctuating  lift 
coefficient  on  the  section,  Sc  is  the  correlation  area  for  the  local  instantaneous 
pressure  fluctuation,  and  U  is  the  relative  mean  velocity  of  the  blade  at  y. 

The  order  of  magnitude  of  the  source  strength  was  suggested  by  Lilley.  It 
was  suggested  that  the  lift  fluctuations  are  related  to  time  fluctuations  in 
the  boundary  layer  thickness  at  the  blade  trailing  edge.  Thus,  it  is  argued 
that  for  typical  fan  conditions,  the  fluctuating  lift  coefficient  should  be  of 
the  order  of  R  -  where  R  is  the  Reynolds  number  based  on  the  cp.ord.  Further 

the  frequency  of  lift  fluctuation  of  the  same  order  as  the  vortex  shedding  fre¬ 
quency,  which  is  approximately  0.2  U/C.  In  addition,  the  correlation  area  which 
is  governed  by  the  large  scale  eddies  is  assumed  to  be  of  the  order  10"^  C/2. 

Using  these  orders  of  magnitude  estimates,  equation  (II-B11)  becomes 


W  =  po 
120tt  c-3 


(II-B12) 


(b)  Noise  from  a  Blade  Operating  in  Turbulent  Flow.  For  a  blade 
operating  in  a  turbulent  flow  a  convenient  form  can  be  obtained  by  introducing 
an  average  lift  curve  slope,  such  that  =  <J>  v/U  (v  being  the  component  of 
turbulent  velocity  normal  to  the  chord) .  Then 


at  u2  at 


Assuming  a  harmonic  form  for  the  turbulent  velocity  component 


where,  to,  is  the  angular  center  frequency,  and  that  S  u>2  is  of  the  order  ot  U1-, 
the  acoustic  power  from  equation  (II-B11)  becomes 
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(XI-B13) 


W  * 
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(c)  Comparison  with  Experimental  Results.  Equations  (II-B12)  and 
(II-B13)  were  developed  for  an  isolated  airfoil.  For  axial  flow  fan  engines 
there  are,  B,  blades  contributing  to  the  overall  acoustic  power,  and  the  blade 
span  is  the  extent  between  the  hub  and  the  tip  of  the  rotor.  Thus  equations 
(II-B12)  and  (II-B13)  may  be  written  as 
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The  total  acoustic  power  levels  calculated  from  the  equations  above  for  a  fan 
with  outlet  guide  vanes  are  shown  in  Figure  (II~B2) .  The  data/theory  comparisons 
are  hared  on  General  Electric  development  vehicle  (D/V  1)  with  the  following 
parameters. 


Number  of  rotor  blades 

38 

Number  of  inlet  guide  vanes 

80 

Chord  of  rotor  blade  (tip) 

4.05  in. 

Chord  of  vane 

1.81  in. 

Average  tip  radius  of  rotor 

17.475  in. 

Average  tip  radius  of  OGV 

It. 975  in. 

Average  hub  radius  of  OGV 

10.325  in. 

Design  tip  speed  of  rotor 

1388  fps 

Average  lift  curve  slope 

0 . 9tt 

The  mean  velocity  relative  to  the  blade  at  each  radial  position  was  taken  as 
the  average  mean  velocity  at  the  leading  edge  and  that  is  the  trailing  edge. 

Hie  measured  and  predicted  values  shown  in  Figure  (IT-B2)  are  in  remarkable 
agreement  with  each  other.  The  sound  power  level  due  to  each  mechanism  is 
piotted  in  Figure  (II-B3),  Vortex  shedding  from  the  rotor  is  seen  to  be  the 
dominant  noise  source. 
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BLADE  TIP  SPEED,  MACH  NO. 


FIGURE  II- B3 


BROADBAND  NOISE  MECHANISMS  ON  A  HIGH  SPEED  FAN 


(3)  Spectral  Distribution  of  Broadband  Noise,  The  text  above  illus¬ 
trated  what  key  aerodynamic  parameters  affect  the  overall  acoustic  sound  power 
and  a  method  for  predicting  the  overall  level  of  noise  due  to  vortex  shedding 
and  due  to  inlet  turbulence.  The  following  will  describe  a  method  which  ex¬ 
tends  the  above  analysis  to  include  the  predicting  spectral  distribution  of 
noise,  and  the  influence  of  aerodynamic  parameters. 

Hence,  consider  the  integral  equation  for  acoustic  pressure  given  by  equation 
(II-B9) 


Noise  resulting  from  Noise  resulting  from  the 

fluctuating  shearing  stresses  effect  of  solid  boundaries 

upon  the  flow  field 

The  first  term  of  the  above  equation  represent  Lighthill's  classical  quadru- 
pole  noise  term,  and  thus  may  be  neglected  since  we  are  considering  dipole  con¬ 
tributions  alone.  In  a  like  vein,  by  dimensional  considerations,  the  first 
two  terms  of  the  sui "ace  integral  may  be  neglected  since  they  represent  monopole 
noise  contributions.  The  resulting  equation  is  then 


P(r,t) 


4"rc 


rr 
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1e 

at 


dS 


(II-B14) 


(a)  Fluid  Dynamics  Related  to  Dipole  Noise.  In  order  to  obtain 
the  relationship  between  the  acoustic  field  and  the  detailed  aerodynamic  para¬ 
meters,  one  can  consider  the  integral  form  of  the  conservation  of  momentum  for 
flow  over  an  isolated  airfoil  as  shown  in  Figure  (II-B4) . 
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(II-B15) 
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FIGURE  II-B4 


CONTROL  VOLUME  FOR  FLOW 
OVER  AN  ISOLATED  AIRFOIL 


where 


V  Control  volume  fixed  in  space 
S  Surface  area  of  control  volume 
n  Unit  outer  normal  vector 
t  Stress  Tensor  (-pi  +  t') 
or  upon  expansion  of  the  stress  tensor 
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Further,  by  assuming  that  the  transverse  components  of  the  fluctuations  of 
mass  flux  are  large  compared  to  the  tangential  components,  equation  (II-B16)  may 
be  written  as 


pu  .  ^  dV  + 


Pl«l 


(-1)  uo 


i  (&2+  u£)  dS2 


i  ds  + 


n)  .  ^  dS 


u,  =  U  i 
1  o- 

u2  =  (U2  +  u^)  i  +  v2  ^ 


where 


Differentiating  with  respect  to  time,  a  form  is  obtained  which  is  directly 
applicable  to  equation  (II-B14)  i.e. 


Dipole  noise  due  Dipole  noise  due  Dipole  noise  due 

to  fluctuations  of  to  fluctuation  of  to  fluctuation  of 

mass  flux  within  momentum  flux  viscous  stresses 

the  control  through  the  wake 

volume  surface  area  S2 
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Hence 


P(r,t) 


_ 1_ 
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(XI-B18) 


Equation  (II-B13)  represents  the  dipole  noise  resulting  from  the  effect  of 
solid  boundaries  upon  the  flow.  The  dipole  noise  thus  represented  is  due  to 
three  phenomenon.  The  first  dipole  source  is  due  to  volumetric  contributions 
to  fluctuations  of  mass  flux,  the  second  dipole  source  is  due  to  the  fluctua¬ 
tions  of  transverse  momentum  flux  through  the  surface  area  of  the  wake,  the 
third  term  is  dipole  noise  due  to  fluctuation  in  the  transverse  components  of 
shear  stress.  In  what  follows,  the  assumption  will  be  made  that  the  shear 
stress  fluctuations  are  negligible  in  comparison  to  the  volume  contributions, 
and  fluctuations  of  momentum  flux. 

(1)  Volume  Contributions  to  Dipole  Noise.  The  total  radiated 
energy  in  the  free  field  for  the  volume  contributions  of  dipole  noise  is,  for 
a  surface  whose  dimensions  are  considered  small  in  comparison  with  wave  length. 


(II-B19) 


Thus  a  correlation  volume  may  be  defined  as 


dV' 


Therefore; 
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or  neglecting  second  order  effects 
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(II-B20) 


where 


Vc  =  Correlation  Volume 

V  =  Wake  Volume 
w 

Further  consider  that  it  is  the  first  term  which  reflects  the  character  of 
the  entire  term  in  the  brackets,  then  equation  (II-B20)  reduces  to 
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or  using  the  isentropic  flow  relations 
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(II-B21) 


where  M  is  the  local  Mach  number  in  the  wake  region,  U  is  the  local  mean  axial 
speed  in  the  wake,  v*  is  the  local  turbulent  transverse  velocity  component  and  y 
is  the  ratio  of  specific  heats. 


The  spectral  distribution  of  acoustic  power  may  be  found  by  taking  the  Fourier 
Transform  of  equation  (IT-B21)  yielding 


W(f)  a,  (p  /a 
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(II-B22) 
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By  assuming  that  the  correlation  volume,  Vc,  is  proportional  to  a  charac¬ 
teristic  dimension  of  the  wake,  to  the  third  power,  namely  the  wake  thickness, 
and  that  the  wake  volume  Vw.  is  proportional  to  Xbfly,  where  X  is  a  reference 
longitudinal  dimension,  b  the  blade  span,  and  Ay  the  wake  thickness,  equation 
(II-B22)  may  be  written  as 


W(f)  * 
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(II-B23) 


Thus,  eauation  (II-B23)  represents  the  spectral  energy  distribution  due  to 
the  volume  contributions  of  the  fluctuations  of  mass  flux.  Further  the  spectral 
distribution  is  seen  to  be  proportional  to  the  fourth  power  of  wake  thickness, 
the  fourth  power  of  frequency,  the  second  power  of  the  transverse  component  of 
turbulent  intensity  and  local  wake  Mach  number,  and  linearly  proportional  to  the 
blade  span. 


(c)  Surface  Contributions  to  Dipole  Noise,  if  0ne  considers  the 
form  of  the  spectral  distribution  of  acoustic  power  due  to  the  surface  contri¬ 
butions  represented  by  the  second  integral  of  equation  (II-B18),  one  can  use 
for  this  limiting  case  the  equation  for  acoustic  pressure 


p(r,t) 
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There  exists  a  correlation  area  S^,  such  that  equation  (II-B24)  may  be 
written  as 
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(II-B25) 


Using  the  expression  for  the  acoustic  power  for  a  dipole,  the  isentropic 
flow  relations,  and  the  Fourier  transform  of  the  resulting  dipole  power  ex¬ 
pression  yields  the  following  spectral  distribution  of  acoustic  energy  due  to 
surface  contributions  for  wake  noise. 
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(II-B26) 


By  assuming  that  the  correlation  area,  Sc  is  proportional  to  the  wake  thickness 
squari  y2,  and  that  the  wake  area,  A„  is  proportional  to  the  span  times  the  wake 
thickness  (bAy)  equation  (II-B26)  becomes 
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(1I-B27) 


Here  it  is  seen  that  the  spectral  acoustic  power  distribution  is  proportional 
tc  wake  thickness  cubed.  Ay-3,  local  Mach  number  to  the  fourth,  M2^,  frequency 
squared,  f2,  and  the  square  of  the  transverse  turbulent  velocity  component  in 
the  wake. 

It  can  be  shown  that  the  acoustic  energy  intensity  and  acoustic  power  for  a 
generalized  dipole  noise  source  varies  as  f^,  i.e.,  for  the  acoustic  power  tt 


n  'v.  4/3  (tt/c)3  f4  pDw 


where  Dy  is  the  strength  of  the  dipole.  Hence  the  formulation  presented  for 
volume  contributions  to  dipole  noise  is  seen  to  check  with  the  frequency  to 
the  fourth  power  law  for  a  radiating  transverse  dipole. 

It  can  be  seen  from  the  above  that  the  airfoil  broadband  noise  gene  "'  '.on 
is  directly  linked  to  the  velocity  fluctuations  around  the  airfoil,  i-  <:  led 
knowledge  of  airfoil  wake  characteristics  is  howevtr  very  restricted.  A  tunnel 
test  program  was  therefore  undertaken  to  permit  a  better  understanding  of  un¬ 
steady  aerodynamic  phenomena  around  airfoils  and  subsequently  broadband  noise 
The  program  is  described  in  ^he  next  section  of  the  report. 
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NOMENCLATURE  SECTION  II-B 
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V,  v' 
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blade  spar 
speed  of  sound 
blade  chord 

local  two-dimensional  lift  coefficient 

frequency 

Green's  Function 

unity  tensor 

unit  vectors  in  longitudinal  and  transverse  directions 

local  Mach  number 

unit  outer  normal  vector 

pressure 

radial  vector  from  origin  or  coordinate  system  to  point  of 
observation  (see  Figure  II-B1) 

position  vector  from  source  location  to  observer  (See  Figure 
II-B1) 

Reynold's  number  based  on  chord  length 

surface  area 

time 

local  longitudinal  mean  velocity  component  in  the  wake 
vector  velocity 

longitudinal  turbulent  velocity  component 
volume 

transverse  turbulent  velocity  component 
acoustic  power 

longitudinal  distance  from  airfoil 
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Y 

A 

6 

X 

U 

it 

P 

'b 

P 

T 

T 

41 

cj 


ratio  of  specific  heats 
increment  of  distance 
Dirac  delta  function 
second  viscosity  coefficient 
first  viscosity  coefficient 
acoustic  power 
density 

acoustic  source  function 
stress  tensor 
viscous  stress  tensor 
average  lift  curve  slope 
angular  frequency 


SUBSCRIPTS 

0  Ambient  conditions,  reference  to  source  coordinate 

system  (See  Figure  II-B1) 

1  Station  one  (See  Figure  II-B4) 

2  Station  two  in  wake  region  (See  Figure  II-B4) 
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(C)  Experimental  Investigation  of  Vortex  Shedding  Noise  in  Subsonic  and 
Transonic  Wind  Tunnels.  Experiments  were  conducted  to  gain  knowledge  of  the 
broadband  noise  generating  mechanisms.  Attention  was  focused  on  vortex  shedding 
noise  which  was  determined  (see  Figure  II-B3)  to  be  the  major  broadband  noise 
contributor.  Basic  experiments  were  conducted  in  subsonic  and  transonic 
wind  tunnels  to  investigate  vortex  bedding  from  simple  geometric  objects  such 
as  cylinders  and  isolated  airfoils.  Data  was  also  obtained  on  loaded  airfoils 
in  two  dimensional  cascades.  Hot  film  anemometers  were  used  to  determine  aero¬ 
dynamic  properties.  Attempts  were  made  to  correlate,  wherever  possible,  the 
velocity  fluctuations  with  the  generated  noise. 

(1)  Test  Facilities. 

(a)  Subsonic  Wind  Tunnel.  The  subsonic  wind  tunnel  is  shown  in 
Figures  J.I-C1  and  II-C2.  The  tunnel  consists  of  a  rectangular,  acoustically 
treated  muffler  and  a  test  section  bolted  on  to  a  29.5  inch  circular  plenum 
chamber.  Interior  dimensions  of  the  muffler  are:  47.5  x  4  x  16  inches  ana  of  the 
test  section  21  x  4  x  16  inches.  (Dimensions  are  given  as;  length,  width  and 
height  respectively) . 

The  plenum  chamber,  acting  as  a  settling  reservoir,  was  provided  with  four 
stainless  steel  screens;  one  of  which  was  a  heavy  gauge  screen  with  about  50% 
open  area.  The  three  others  contained  fine  20  mesh  screening.  The  muffler  was 
incorporated  in  the  system  to  eliminate  upstream  piping  noise.  The  acoustic 
treatment  linings  on  two  sidewalls  of  the  muffler  section  consisted  of  1  inch 
thick  Scottfelt,  located  between  the  muffler  wall  and  a  stainless  steel  perforated 
plate  (22%  open  area,  0.03  inches  thick).  Before  the  flow  enters  the  test,  section 
it  meets  two  contraction  blocks  which  form  a  bellmouth  inlet  to  the  test  section. 
The  air  is  exhausted  to  atmospheric  conditions.  The  exhaust  cross  sectional 
area  was  larger  than  the  cross  sectional  area  of  the  test  section  (5/1  expansion 
area)  so  that  the  exhausted  air  could  expand. 

(b)  Transonic  Wind  T1  rnel .  A  cross  section  of  the  transonic  wind 
tunnel  facility  is  shown  in  Figure  IT.-Cu.  Inside  the  tank  a  cascade,  placed 
between  two  Plexiglas  sidewalls,  was  mounted  in  a  moveable  drum.  This  drum 

can  rotate  in  relation  to  the  oncoming  flow,  thereby  varying  the  inlet  air  angle 
to  the  cascade.  The  tunnel  was  ea.uipped  with  instrumentation  to  obtain  aerodynamic 
and  turbulence  data.  Upstream  and  downstream  pressure  and  hot  film  traverses 
were  executed  using  actuators.  The  actuators  are  remotely  controlled,  and  are 
movable  in  all  directions.  Noise  data  was  recorded  on  a  microphone  located 
6  feet  away  from  the  cascade  inside  the  wind  tunnel. 

A  bloc’  diagram  of  the  piping  system  for  the  transonic  wind  tunnel  is  shown  in 
Figure  II-C4.  The  test  section  of  the  wind  tunnel  is  at  the  end  of  the  pipe 
entering  the  test  tank.  Atmospheric  air  enters  the  system  at  (A)  and  can  be  routed 
through  the  air  dryer  or  sent  directly  to  the  test  tank.  This  inlet  air  then  passes 
through  a  plenum  chamber  and  a  bellmouth  before  entering  the  test  section.  After 
exhausting  into  the  tank  the  air  is  drawn  out  the  bottom  of  the  tank  through 
valve  (C)  to  the  inlet  of  the  compressor.  Valve  (C)  regulates  the  flow  velocity 
by  controlling  the  back  pressure  in  the  tank.-  From  the  .ischarge  of  the  com¬ 
pressor  most  of  the  air  :s  exhausted  to  the  atmosphere  but  a  portion  of  it  is 
recirculated  through  the  system  by  valve  (B) . 
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rIGURE  II-C1.  MAJOR  COMPONENTS  OF  THE  SUBSONIC  WIND  TUNNEL 


FIGURE  II- C2  PHOTOGRAPH  OF  THE  SUBSONIC  WIND  TUNNEL 


FIGURE  II-C3 .  CROSS-SECTION  01  TRANSONIC  CASCADE  TUNNEL 


ATMOSPHERE 


FIGURE  II-C4.  BLOCK  DIAGRAM  OF  THE  TRANSONIC  WIND  TUNNEL 


Velocities  in  the  tunnel  are  set  by  using  a  mechanical  pressure  gauge  to 
record  plenum  chamber  total  pressure  and  static  pressure  upstream  of  the  test 
section  and  forming  the  ratio  Ps/P0  where  P$  is  static  pressure  and  P0  is 
total  pressure.  Manometer  banks  filled  with  TBE  manometer  fluid  are  provided 
to  measure  an  upstream  and  downstream  pressure. 

(2)  Instrumentation. 

(a)  Hot  Film  Anemometers  -  Velocity  Measurements.  It  is  advantageous 
to  obtain  velocity  data  by  hot  film  probes  since  the  compactness  of  such  probes 
minimizes  the  disturbance  to  the  flow  and  also  because  of  the  fast  response 

time  of  the  probe  to  flow  changes.  The  hot  film  sensors  are  an  improved  version 
of  the  hot  wire  in  that  they  are  more  rigidly  constructed.  The  hot  film 
sensor  consists  of  a  thin  film  of  metal  deposited  onto  a  ceramic  substrate  held 
by  needle  shaped  supports  which  in  turn  are  attached  to  the  stem  of  the  probe. 

An  X-array  of  sensors  and  probe  assembly  is  shown  in  Figure  II-C5.  Sensors, 
probes  and  associated  instrumentation  were  purchased  from  Thermo  Systems,  Inc. 

In  operation,  each  sensor  is  part  of  a  Wheatstone  bridge.  The  sensors  are 
kept  at  constant  temperature  by  the  aneomometers ,  shown  in  the  block  diagram  of 
Figure  II-C6.  The  equipment  used  is  shown  on  Figure  II-C7.  They  instantaneously 
measure  fluid  flow  parameters,  by  sensing  the  neat  transfer  rate  between  the 
electrically  heated  sensor  and  the  flow.  When  the  flow  velocity  changes,  com¬ 
pensation  is  provided  to  maintain  a  constant  temperature  difference.  Since  the 
resistances  in  the  Wheatstone  bridge  are  kept  constant,  through  a  feedback  control 
in  the  circuitry,  the  voltage  across  the  bridge  is  directly  proportional  to  the 
current,  I,  through  the  sensor  and  the  power,  P,  is  equal  to  the  product  of  the 
square  of  the  current  and  the  resistance,  R  (P  =  I2  R).  Therefore,  the  square 
of  the  voltage,  measured  on  top  of  the  bridge  is  directly  proportional  to  the 
instantaneous  heat  transfer  between  the  sensor  and  its  environment.  This 
heat  transfer  is  proportional  to  the  fluid  velocity. 

The  output  of  the  sensor  amplified  by  the  bridge  is  not  linear  with  flow 
velocity  changes;  therefore,  a  linearizer  was  incorporated  in  the  system  (See 
Figure  II-C6  and  II-C7)  for  convenience.  An  example  of  a  linearized  calibration 
curve  for  an  X-array  sensor  is  shown  in  Figure  II-C8.  To  obtain  calibration 
of  sensors  and  flow  parameters  a  correlator  was  used.  This  device  is  able 
to  perform  additions  and  subtractions  or  alternatively  present  the  output  of  the 
sensors  to  a  DC  Voltmeter  and  a  set  of  four  "true  root  mean  square"  voltmeters 
(TRI1S) .  These  voltmeters  display  mean,  mean  square,  and  RMS  values  of  the  sensor 
outputs.  The  voltmeters  were  altered  to  merge  with  a  Honeywell  digitizing  system. 

This  system  consists  of  a  multipoint  analog  scanner  which  sequentially  scans  the 
outputs  of  the  various  voltmeters.  These  outputs  are  then  sequentially  processed, 
displayed  on  the  digital  voltmeter,  converted  to  digital  form  and  fed  to  the  teletype. 
It  was  convenient  to  do  so,  since  the  teletype  provided  a  printed  copy  of  the  data 
and  a  punched  paper  tape. 

(b)  Comk  user  Microphone  -  Sound  Measurements .  Sound  measurements 
were  obtained  during  testing  in  the  cascade  using  a  1/4  inch  Bruel  and  Kjaer  (B&K) 
condenser  microphone.  The  microphone  was  located  6  feet  away  from  the  exit 
measuring  station  in  the  cascade  tunnel.  Long  extension  cables  made  it  neces¬ 
sary  to  use’  a  power  supply  (BhK  2801)  in  conjunction  with  the  microphone  (See  Figure 
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FIGURE  II-C5.  VELOCITY  SENSITIVE  HOT  FILM  PROBE 
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FIGURE  II-C6.  BLOCK  DIAGRAM  OF  THE  HOT  FILM  INSTRUMENTATION 
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FIGURE  II-C7  PHOTOGRAPH  OF  THE  HOT  FILM  INSTRUMENTATION 


OUTPUT  OF  SENSOR  IN  VOLTS 


FIGURE  II-C8.  CALIBRATION  CURVE:  OUTPUT  OF 
SENSOR  VS.  FLOW  MACH  NUMBER 


II-C9).  This  apparatus  supplied  ample  voltage  to  the  microphone's  cathode 
follower  and  also  provided  the  necessary  polarization  voltage  to  the  condenser 
microphone.  Calibration  of  the  microphone  was  accomplished  by  temporarily 
attaching  a  B&K  Pistonphone  (Type  4220)  to  the  microphone.  The  sound  pressure 
level  of  this  source  was  124  dB  within  ±0.?  dB. 

The  output  signal  of  the  microphone  was  connected  to  a  B&K  Wave  Analyzer 
(Type  2107).  This  analyzer  utilizes  a  6%  constant  percentage  bandwidth.  It 
performed  on  line  to  produce  narrowband  analysis  of  the  noise.  A  built-in 
mechanical  device  (combination  of  gears)  provided  automatic  turning  when 
coupled  to  the  motor  of  the  B&K  level  Recorder  (Type  2305) .  Combined 
in  this  manner,  frequency  amplitude,  diagrams  of  the  sound  spectrum  were  obtained. 

An  alternate  routing  was  made  using  a  Lockheed  411B,  4  channel  AM  Tape  Recorder 
which  preserved  the  noise  signal  for  off  line  narrowband  analysis. 

(3)  Subsonic  Wind  Tunnel  Tests.  The  inlet  turbulence  and  the  inlet 
velocity  profile  were  checked  at  M  =  0.4  using  a  hot  film  traverse.  The 
traverse  originated  at  the  center  of  the  duct,  represented  by  Y  =  0,  and 
ended  near  the  top  wall,  indicated  by  Y  =  8.  Data  reduction  of  the  hot  film 
output  yi£lded  turbulence  intensity  In  the  axial  flow  direction  (u'/U)  and  the  mean 
velocity  U  for  each  probe  position.  Figure  1I-C10  shows  the  mean  velocity  profile 
for  the  air  flowing  through  the  duct.  The  shape  of  the  mean  velocity  profile 
was  flat  except  at  Y  =  5  inches  where  the  mean  velocity  decreased  due  to  the 
influence  of  the  solid  boundary  (top  wall).  As  can  be  seen  in  Figure  II-C11 
the  turbulence  intensity  u'/U  in  the  center  of  the  duct  is  less  than  4%.  About 
one  inch  away  from  the  top  wall  the  u'/U  =  7.7%. 

(a)  Cylinder  Wake  Investigation.  The  wake  of  a  one-inch  cylinder 
was  explored  by  means  of  the  hot  film  anemometers,  described  in  Section  Il-C2a. 
Turbulence  intensity  measurements  were  made  for  upstream  and  downstream  stations 
relative  to  the  cylinder  with  flow  Mach  numbers  0.1,  0.2  and  0.3. 

Figure  II-C12  depicts  the  cest  section  in  whidi  the  cylinder  was  mounted. 

Figure  II-C13  shows  the  turbulence  intensity  u'/U  on  the  wake  axis  of  the 
cylinder  as  a  function  of  downstream  distence  for  various  Reynolds  numbers. 

As  expected,  the  turbulence  intensity  decreased  for  increasing  downstream  distances. 
The  turbulence  intensity  increased  for  increasing  flow  Mach  number.  While  Figure 
H-C13  indicates  the  turbulence  intensity  normalized  to  the  local  mean  velocity 
(U)  ;  Figure  II-C14  shows  the  turbulence  intensity  normalized  to  the  undisturbed 
(i.e.,  upstream)  mean  velocity  (U  )  of  the  flow  versus  downstream  positions.  The 
results  were  compared  to  experiments  of  Kovasnay  (Reference  1)  and  show  that  there 
is  a  point  on  the  wake  axis  of  the  cylinder  where  the  turbulence  intensity  peaks, 
results  indicate  this  point  to  be  5  diameters  downstream  of  the  cylinder  while 
lower  Reynold's  number  experiments  put  in  at  7. 

Figure  II-C15  shows  the  turbulence  intensities  in  the  flow  direction  (u'/U) 
two  inches  upstream  of  the  cylinder  for  several  positions  in  the  duct  and  for 
three  flow  Mach  numbers.  The  turbulence  intensity  is  Jess  than  5%  in  the  center 
of  the  duct  and  is  shown  to  increase  rapidly  when  approaching  the  duct  walls. 


FIGURE  II-C9.  DATA  ACQUISITION  BLOCK  DIAGRAM. 


76 


IMMERSION  Y,  INCHES 


FIGURE  II-C11.  INLET  TURBULENCE  OF  THE 

SUBSONIC  WIND  TUNNEL. 
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FIGURE  II-C12.  TEST  SECTION  WITH  UPSTREAM  AND 
DOWNSTREAM  MEASURING  STATIONS. 
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FIGURE  II-C13.  TURBULENCE  INTENSITY,  NORMALIZED 
TO  LOCAL  MEAN  VELOCITY,  OF  THE 
WAKE  BEHIND  A  CYLINDER. 
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FIGURE  II-C14.  TURBULENCE  INTENSITY,  NORMALIZED  TO  UNDISTURBED 

(UPSTREAM)  MEAN  VELOCITY,  OF  THE  WAKE  BEHIND 
A  CYLINDER. 
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FIGURE  II-C15.  UPSTREAM  TURBULENCE  INTENSITY 
VS.  IMMERSION. 
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In  addition  to  turbulence,  noise  and  vortex  shedding  frequency  determinations 
were  made  for  a  1/2  inch  cylinder  for  three  upstream  flow  Mach  numbers.  Figures 
II-C16,  18  and  19  give  the  turbulence  intensities  in  the  direction  of  the  flow 
for  three  downstream  stations  versus  probe  position  in  the  duct  for  three  flow  Mach 
numbers.  Figure  II-C20,  21  and  22  show  the  turbulence  intensities  w'/U  (in 
the  spanwise  direction  of  the  cylinder  )  at  three  downstream  stations  for 
three  flow  Mach  numbers.  All  these  figures  show  that  the  turbulence  intensities 
decrease  with  increasing  flow  Mach  number.  This  can  be  explained  by  considering 
the  hot  film  probe  to  the  positioned  at  one  point  in  the  wake  of  the  cylinder. 

If  the  flow  Mach  number  is  increased, the  output  of  the  hot  film  increases,  due  to 
the  linear  response  of  the  hot  film  with  flow.  Furthermore;  the  mean  velocity  of 
the  flow  corresponds  to  the  mean  value  of  the  output  of  the  hot  film,  while 
the  turbulence  u*  and  w*  corresponds  to  small  signal  fluctuations  about  this 
mean  le-'el.  Hence,  when  the  flow  Mach  number  increases  the  mean  value  of  the 
hot  film  output  increases.  Since  this  value  appears  in  the  denominator  of  the 
turbulence  intensity  it  is  evident  that  the  turbulence  intensity  is  decreasing 
for  increasing  flow  Mach  numbers.  Figure  1I-C17  shows  the  calculated  turbulence 
u'  on  the  wake  axis  of  the  cylinder,  using  the  turbulence  intensities  shown 
in  Figure  II-C16.  As  can  be  seen  the  u*  turbulence  increased  with  increasing  flow 
Mach  number. 

Using  a  single  hot  film  a  radial  traverse  was  made  of  the  half  inch  cylinder 
wake  and  the  probe  stopped  at  the  position  of  maximum  turbulence  (position 
where  vortices  are  shed  off  the  cylinder) .  This  output  was  fed  into  a  General 
Radio  wave  analyzer  for  spectral  analysis.  The  frequency  where  the  signal 
was  a  maximum  was  noted  to  be  the  vortex  shedding  frequency.  Figures  II-C23, 

24,  25  show  the  frequency  spectra  for  the  half-inch  cylinder  for  flow  Mach 
numbers  0.1,  0.2  and  0.3  respectively.  The  position  for  which  the  turbulence 
was  a  maximum  was  found  to  be  0.25  inches  from  the  cylinder's  axis  at  two 
diameters  downstream  of  the  cylinder. 

(b)  Noise  Measurements  on  a  Cylinder.  The  sound  radiated  by  a 
cylinder  immersed  in  an  airstream  has  been  studied  by  a  number  of  investigators 
(Reference  2).  Etkin  (Reference  3)  found  that  the  pure  tone  in  the  sound  spectrum 
was  associated  with  the  shedding  of  vortices  .rom  the  cylinder.  The  frequencies 
were  found  to  be  the  same. 

In  considering  the  sound  intensity  spectrum,  investigators  found  that: 

The  pure  tone  associated  with  the  vortex  shedding  frequency 
of  the  cylinder  is  related  to  the  alternating  lift  force  acting 
on  the  cylinder  and  has  the  directionality  of  a  line  dipole  with 
the  axis  normal  to  the  flow. 

The  tone  generated  at  twice  the  vortex  shedding  frequency  is 
related  to  the  fluctuating  drag,  with  the  directionality  of  a 
line  dipole  with  axis  parallel  to  tiu  stream. 

The  acoustic  power  radiation  varies  as  u”,  where  n  lies  between  5 
and  6  depending  on  the  ratio  of  cylinder  length  to  cylinder  dia¬ 
meter,  i.e.,  £/d. 
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FIGURE  II-C16,  DOWNSTREAM  TURBULENCE  INTENSITY  VS  IMMERSION, 
TWO  DIAMETERS  DOWNSTREAM  OF  A  CYLINDER, 
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URBULENCE  INTENSITY  IN  PERCENT 


FIGURE  II-C18  „  TURBULENCE  INTENSITY  VS  IMMERSION,  TEN 
DIAMETERS  DOWNSTREAM  OF  A  CYLINDER. 
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FIGURE  II-CI9.  TURBULENCE  INTENSITY  VS  IMMERSION, 
EIGHTEEN  DIAMETERS  DOWNSTREAM  OF  A 
CYLINDER. 
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FIGURE  II-C20.  TURBULENCE  INTENSITY  VS  IMMERSION, 

TITO  DIAMETERS  DOWNSTREAM  OF  A  CYLINDER 
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FIGURE  II-C22.  TURBULENCE  INTENSITY  VS  IMMERSION 
EIGHTEEN  DIAMETERS  DOWNSTREAM  OF 
A  CYLINDER. 
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FIGURE  II-C25.  VORTEX  SHEDDING  FREQUENCY  FOR  1/2”  CYLINDER 

TWO  DIAMETERS  DOWNSTREAM. 
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The  main  limitations  o£  the  theory  are  that  one  must  know  the  forces  acting 
on  the  cylinder  and  the  mechanism  for  converting  aerodynamic  to  acoustic  energy. 

An  acoustic  probe  system  was  placed  inside  the  wind  tunnel  at  the  end  of  the 
test  section,  i.e.,  downstream  of  the  cylinder  in  a  plane  that  bisects  the 
cylinder  span.  The  system  was  calibrated  and  the  probe  system  losses  determined 
in  each  frequency  band.  The  microphone  signal  was  fed  directly  into  the  B&K 
analyzer . 

In  addition  two  microphones  were  employed:  a  1/2  inch  microphone  16"  down¬ 
stream  of  the  cylinder  and  another  microphone  (1/4  inch)  8"  /artically  below  the 
cylinder.  These  two  microphones  were  flush  mounted  with  the  wind  tunnel  walls, 
i.e.,  only  the  microphone  diaphragms  were  exposed  to  the  air  stream. 

Figures  II-C26  through  II-C37  show  the  sound  spectra  (SPL  versus  frequency, 

"f")  at  various  flow  Mach  numbers,  of  the  noise  with  and  without  the  cylinder  at 
the  two  microphone  locations.  The  SPL  differences  between  the  background  noise 
(noise  without  the  cylinder)  and  noise  with  the  cylinder  in  place  are  given  in 
Figure  II-C28,  31,  34  and  37.  In  these  figures  the  frequency  of  the  maximum  dif¬ 
ference  between  noise  with  and  without  the  cylinder  is  clearly  visible  and  indicated 
by  ft. 

Table  II-C1  shows  the  measured  vortex  shedding  frequency  fv,  as  determined 
with  the  hot  film  probe  in  the  wake  of  1/2  inch  cylinder,  and  the  tone  frequency 
ft  determined  from  the  noise  spectra  measured  by  the  microphones  for  the  three 
Mach  numbers.  The  Strouhal  number  (Sh)  was  calculated  from  the  data  according 
to  the  relation: 


Sh  =  f  d/U 
v  o 

It  can  be  seen  that  the  vortex  shedding  frequency  as  measured  by  the  hot  film 
and  the  noise  peak  as  measured  by  the  microphones  are  nearly  the  same.  Thus, 
the  maximum  turbulence  is  related  to  the  maximum  noise. 

Figure  II-C38  shows  the  OASPL  versus  flow  Mach  number  for  a  half  inch  cylinder. 
The  CSArL  was  read  on  the  General  Radio  Sound  Level  Meter.  The  data  represented 
by  the  short  dashed  curve  was  obtained  from  the  microphone  located  underneath 
the  cylinder  while  the  solid  curve  was  obtained  from  the  microphone  downstream 
of  the  cylinder.  The  slope  of  these  curves  was  equal  to  2.5  for  the  microphone 
below  the  cylinder  and  3.0  for  the  microphone  downstream  of  the  cylinder.  Figure 
II-C39  shows  the  SPL  at  the  vortex  shedding  frequency  generated  by  the  cylinder 
for  various  flow  Mach  numbers. 

Making  use  of  the  following  relationships: 
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ONE-HALF  INCH  CYLINDER 

M 

£  :  Hot  Film 

V  (Hz) 

ft:  1/2"  Mike 

Below  Cylinder 
(Hz) 

Sh 

0.1 

520 

500 

0.19 

0.2 

990 

1000 

0.18 

0.3 

1480 

1250 

0.18 

0.4 

X 

2000 

X 

X  =  Not  Available 


TABLE  II-C1 


Comparison  Between  Hot  Film  and  Microphone  Vortex  Shedding  Frequencies 


FIGURE  II-C26.  1/3  OCTAVE  SOUND  SPECTRUM  OF  A  1/2"  CYLINDER  IMMERSED 

IN  AN  AIR  STREAM. 
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FIGURE  XI-C27.  1/3  OCTAVE  SOUND  SPEC'  ROM  OF  A  1/2"  CYLINDER  IMMERSED 
IN  AN  AIR  STREAM. 


140 


100 


-C30.  1/3  OCTAVE  SOUND  SPECTRUM  OF  A  1/2"  CYLINDER  IMMERSED  IN 
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FIGURE  II-C32.  1/3  OCTAVE  SOUND  SPECTRUM  OF  A 

AN  AIR  STREAM. 
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FIGURE  II-C37 .  1/3  OCTAVE  SOUND  SPECTRUM  OF  A  1/2"  CYLINDER  IMMERSED  IN  AN  AIR  S 
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FIGURE  II-C38 .  OASPL  OF  ] 

VS  FLOW  Mi 


and 


SPL  =  10  log  I/Iref  +  0.2, 


it  can  be  shown  that  the  slope  of  the  undisturbed  velocity  curve  can  be 
expressed  as : 


SPL1  -  SPL2 
10  log  Vj/V” 


where  "a"  is  a  constant,  "U"  the  undisturbed  velocity  of  air,  and  I  the  intensity 
of  sound.  Etkin's  (Ref.  3;  shows  that  the  value  of  n  is  between  5  and  6.  In  the 
present  experiments  the  slope  of  the  curves  is  somewhat  greater.  Below  M= 

0.2  the  slope  was  equal  to  1.6,  this  is  due  to  the  low  signal  to  noise  ratio 
in  the  subsonic  wind  tunnel  at  this  low  flow  Mach  numbers. 

(c)  Symmetric  Airfoil  Wake  Investigation.  The  wake  of  a 
symmetric  airfoil  (inch  chord,  4  inch  span,  25%  thick)  at  zero  angle  of  attack 
was  explored  at  four  flow  Mach  numbers  of  0.1,  \2,  0.3  and  0,4.  Velocity 
and  turbulence  intensity  measurements  were  made  at  both  upstream  and  downstream 
stations.  Figure  II-C40  depicts  the  test  section  in  which  the  airfoil  was 
mounted  along  with  the  distances  between  the  various  measuring  stations.  During 
tests,  the  turbulence  intensity  in  the  axial  or  X  -  direction,  u'/U  was  mea¬ 
sured,  together  with  the  turbulence  intensity  in  the  spanwise  or  Z  -  direction, 
w'/U.  In  Figure  II-C40  the  direction  of  the  fluctuating  velocity  components 
are  drawn  at  point  A.  The  lot  al  mean  velocity  is  denoted  by  U  . 

Figure  II-C41  presents  an  overall  picture  of  the  turbulence  intensity  u'/U 
on  the  wake  axis  as  a  function  of  downstream  distance,  expressed  as  the 
dimensionless  quantity  x/c,  (ratio  of  downstream  distance  to  airfoil  chord). 

Far  away  from  the  airfoil  turbulence  intensities  comparable  to  the  inlet  turbulence 
intensities  are  reached.  As  can  be  seen,  the  turbulence  intensities  decrease 
for  increasing  Mach  number;  however,  the  u'  turbulence  values  are  increasing. 

The  u'  value  was  obtained  by  averaging  the  output  of  the  sensor  for  3  seconds. 

In  Figure  II-C42  the  u'  turbulence  on  the  wake  axis  two  chord  lengths  downstream 
of  the  airfoil  versus  four  flow  Mach  numbers  is  plotted  and  is  seen  to  be 
increasing. 

Figure  II-C43  and  44  are  plots  of  the  turbulence  intensity  u'/U,  at  two 
downstream  stations,  for  four  flow  Mach  numbers  (as  a  function  of  immersion  y/c) . 

As  may  be  noted,  the  intensities  decrease  rapidly  beyond  a  distance  of  3 
airfoil  chords.  Approximately  11  inches  downstream,  the  effect  of  the  airfoil 
on  turbulence  is  hardly  evident.  Thus,  the  axial  extent  of  the  wake  of  this 
airfoil  is,  for  all  practical  purposes,  3  airfoil  chords.  This  information  is 
important  in  establishing  the  correlation  volume  in  the  analysis  of  broadband 
noise  generation  (See  Section  II-B2) . 
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FIGURE  II-C40.  TEST  SECTION  WITH  UPSTREAM  AND 
DOWNSTREAM  MEASURING  STATIONS. 
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FIGURE  II-C4I.  TURBULENCE  INTENSITY  ON  THE  WAKE  AXIS  OF  i 

SYMMETRIC  AIRFOIL  VS.  DOWNSTREAM  DISTANCE 
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FLOW  MACH  NUMBER 


FIGURE  II-C42.  TURBULENCE  ON  THE  WAKE  AXIS,  TWO  CHORD 

LENGTHS  DOWNSTREAM  OF  A  SYMMETRIC  AIRFOIL 
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FIGURE  II-C43.  TURBULENCE  INTENSITY  OF  A 

SYMMETRIC  AIRFOIL  VS  IMMERSION 
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In  Figures  II-C45  and  46  the  turbulence^  intensity  w’/U  is  plotted  as  a 
function  of  immersion  y/c.  Comparing  u’/U  at  the  same  station,  for  example. 

Figures  II-C43  and  46  it  is  evident  that  w*  is  smaller  than  u’. 

An  overall  picture  of  w'/U  versus  downstream  distance  is  given  in  Figure 
III-C47. 

Thus,  in  the  two  dimensional  case  being  examined  the  u*  turbulence  is  the 
dominant  turbulence  component,  c-cd,  as  in  the  case  of  the  cylinder  the  wake 
turbulence  intensity  increases  with  Mach  number  relative  to  the  airfoil. 

(4)  Transonic  Wind  Tunnel  Tests  of  Thin  Airfoils.  Tests  were  carried 
out  with  a  cascade  of  thin  airfoils  (chord  length:  2.7  inches,  t/c:  2.2%) 
in  a  transonic  wind  tunnel  facility.  A  sketch  of  this  tunnel  with  its  assembled 
blade  package  is  shown  in  Figure  II-C3. 

The  incoming  air  velocity  was  varied  from  452  ft/sec  to  891  ft/sec  while 
the  turbulence  intensities  were  measured  one  chord  length  downstream  of  the  cas¬ 
cade.  Turbulence  measurements  in  both  the  mean  flow  direction,  u'/U  and  in  the 
direction  perpendicular  to  the  span  of  the  airfoils  v'/U  were  recorded.  As 
in  t  previous  subsonic  isolated  airfoil  tests,  hot  film  sensors  were  utilized. 

In  audition  sound  measurements  were  made  during  testing  in  the  cascade  using 
a  1/2  inch  Bruel  and  Kjaer  condenser  microphone.  The  microphone  was  located 
6  feet  from  the  exit  measuring  station  in  the  cascade  tunnel. 

Before  examining  the  wake  of  the  airfoils  the  magnitude  of  the  inlet 
turbulence  intensities  i.e.,  the  intensities  one  chord  length  upstream  of  the 
cascade  were  determined.  The  inlet  turbulence  intensities  were  recorded  in  the 
center  of  the  four  inch  wide  rectangular  cascade  duct  while  varying  the  inlet 
Mach  number  (M) .  A f  M  =  0.4,  0.5,  0.7  the  corresponding  u'/U's  were  1,3,  1.1, 

0.7  percent  and  v'/U's  1.2,  1.1,  0.7  percent  respectively.  These  values  are 
small  enough  so  that  the  turbulence  intensities  in  the  wakes  of  the  airfoils 
can  be  distinguished  from  the  tunnel  turbulence  levels.  A  total  pressure  profile 
measurement  was  also  taken  across  the  width  of  the  duct,  upstream  of  the  cascade 
at  M  =  0.7,  using  a  Pitot  tube.  Result;  are  plotted  as  the  ratio  of  the  total 
plenum  pressure  to  the  inlet  total  pressure  in  Figure  II-C48.  The  pressure  profile 
is  flat,  away  from  the  walls. 

The  cascade  consisted  of  11  airfoils  capable  of  being  placed  at  various  angles 
of  attack,  corresponding  to  various  incidence  angles.  The  incidence  angle  beirg 
defined  as  the  angle  between  the  tangent  to  the  airfoil's  mean  camber  line  at-  the 
leading  edge,  and  the  inlet  air  direction.  Characteristics  of  this  cascade  are 
shown  in  Figure  II-C49. 

Figure  I1-C50,  51,  and  52  show  the  mean  axial  velocity  distribution  behind 
the  cascade  for,  respectively,  free  stream  Mach  numbers  of  .4,  ,5,  and  .7.  As 
would  be  expected  the  velocity  defects  increase  with  increasing  Mach  numbers. 

Or  what  is  more  important  from  the  acoustic  point  of  view,  the  wake  volume 
increases  (see  Section  1I-B2). 
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TURBULENCE  INTENSITY  w'/H  IN  PERCENT 


FIGURE  II-C46.  TURBULENCE  INTENSITY  OF  A 

SYMMETRIC  AIRFOIL  VS.  IMMERSION 
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SYMMETRIC  AIRFOIL 
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FIGURE  II-C48 .  TOTAL  PRESSURE  PROFILE  ACROSS  THE  DUCT 
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FIGURE  II-C50.  MEAN  VELOCITY  PROFILE,  ONE  CHORD  DOWNSTREAM 

OF  CASCADE . 


FIGURE  II-C51.  MEAN  VELOCITY  PROFILE,  ONE  CHORD  DOWNSTREAM  OF  CASCADE 


FIGURE  II-C52 .  MEAN  VELOCITY  PROFILE,  ONE  CHORD  DOWNS 
OF  CASCADE. 


The  next  6  figures  show  the  velocity  turbulences.  As  has  been  noted 
previously,  the  turbulence  intensity  decreases  with  increasing  Mach  number; 
however,  the  absolute  turbulence  level  increases.  Figures  II-C59  through  62 
show  che  turbulence  levels  at  .4  and  .7  Mach  numbers.  The  increase  in  u'  and 
v'  turbulences  with  Mach  number  is  clearly  seen. 

Figure  II-C63  is  an  overlay  of  the  noise  data  at  .4  and  .7  Mach  number. 

The  broadband  noise  has  increased  over  most  of  the  audible  frequency  range  at 
the  higher  Mach  number.  Figures  II-C64  and  65  contain  the  corresponding 
turbulence  components  u'  and  v'  at  both  Mach  numbers.  The  increase  in  both 
turbulence  components  with  Mach  number  can  be  seen. 

It  is  interesting  to  note  that  at  12  KHz,  there  is  a  peak  in  both  the 
turbulence  (Figures  II-C64  and  65)  and  noise  (Figure  II-C63)  spectra.  This 
peak  is  believed  to  be  associated  with  the  fluctuation  of  the  airfoil's  lift 
vector. 

Another  parameter  which  is  important  to  broadband  noise  generation  is 
incidence  angle.  Figure  II-C66  shows  the  mean  axial  velocity  distribution  at 
a  Mach  number  of  .7  at  two  different  incidence  angles.  It  is  obvious  that  the 
increased  incidence  angle  has  deepened  the  wake  velocity  defect.  Figures  II-C67 
and  68  show  the  turbulence  spectra  at  3  incidence  angles.  In  general,  the 
turbulence  has  increased  all  across  the  spectrum  with  increasing  incidence  angle. 

Figure  II-C69  is  an  overlay  of  the  tunnel  noise  (no  cascade),  noise  at  an 
incidence  angle  of  5°,  and  noise  at  an  incidence  angle  of  7.5°.  The  background 
noise  due  to  the  facility  is  at  least  10  dB  below  the  noise  with  the  cascade 
package  in  place.  The  7.5°  incidence  angle  noise  is  generally  above  the  5° 
incidence  angle  data;  this  is  particularly  true  at  12  KHz.  Looking  back  to  the 
turbulence  data  of  Figures  II-C67  and  68,  it  can  be  seen  that  a  large  turbulence 
increase  also  occurred  at  12  KHz.  This  is  due  to  the  increased  lift  being 
generated  at  the  higher  incidence  angle  which,  in  turn,  results  in  higher  noise 
and  turbulence  at  the  characteristic  lift  fluctuation  frequency. 

(5)  Summary  of  Wind  Tunnel  Results.  In  both  subsonic  isolated  element 
tests  and  in  cascade  tests,  the  increased  generation  of  broadband  noise  with 
increasing  relative  Mach  number  was  clearly  demonstrated.  In  addition,  the 
content  of  the  noise  spectra  and  the  wake  turbulence  spectra  were  shown  to  be 
related  in  so  far  as  peaks  in  the  turbulence  spectra  occurred  at  the  same 
frequency  as  peaks  in  the  noise  spectra. 

Utilizing  a  cascade  of  airfoils,  it  was  shown  that  the  noise  and  turbulence 
levels  increased  with  increasing  incidence  angle. 

These  results  suggest  that  two  ways  of  reducing  broadband  noise  are  to 
reduce  relative  Mach  numbers  and/or  bring  incidence  angles  closer  to  optimum 
(minimum  wake  turbulence)  angles. 
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(6)  Rotating  Wake  Measurements.  In  order  to  better  understand  the 
viscous  wake  properties  which  affect  noise  generation  on  both  rotating  and 
non- rotating  blades,  the  hot  film  instrumentation  previously  used  in  the  cascade 
was  mounted  on  a  General  Electric  high  bypass  ratio  single  stage  fan  development 
vehicle.  Traverses  were  made  behind  the  rotor  about  one  rotor  chord  downstream 
and  behind  the  OGV's.  The  immersion  positions  are  shown  in  Table  II-C2. 


Position 

Number 


Radial  Immersion 
Depth,  Inches 


1 

2 

3 

4 

5 


.53 

1.64 

2.83 

4.12 

5.56 


TABLE  II-C2 

Probe  Positions  and  Immersion  Depths 
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FIGURE  II-C55.  TURBULENCE  INTENSITY  VS  POSITION  IN  CASCADE  TRANSONIC 
FOR  A  TIP  SECTION  AIRFOIL  ON  D/V  #4 


FIGURE  II-C56.  TURBULENCE  INTENSITY  VS  FO.  VIION  IN  CASCADE  TRANSONIC  TUNNEL 

FOR  A  TIP  SECTION  AIRFOIL  ON  D/V  #4 
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FIGURE  II-C58.  TURBULENCE  INTENSITY  VS  POSITION  IN  CASCADE  TRANSONIC  TUNNEL 

FOR  A  TIP  SECTION  AIRFOIL  ON  D/V  #4 
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FIGURE  II-C60.  u'  TURBULENCE  INTENSITY  DOWNSTREAM  OF  CASCADE. 

FOR  A  TIP  SECTION  AIRFOIL  ON  D/V  #4 
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FIGURE  II-C61.  v'  TURBULENCE  INTENSITY  DOWNSTREAM  OF  CASCADE 

FOR  A  TIP  SECTION  AIRFOIL  ON  D/V  #4 
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FIGURE  II-C65.  v'  TURBULENCE  SPECTRA  FOR  A  TIP  SECTION  AIRFOIL  OF  D/V  #4 
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FIGURE  II-C67.  NARROWBAND  FREQUENCY  ANALYSIS  OF  THE  TURBULENCE  COMPONENT  u’  FOR  A  TIP  SECTION  AIRFOIL  OF  D/V  #4 
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Data  reduction  concentrated  on  narrowband  (50  Hz  bandwidth)  analysis  of 
axial,  u’,  and  tangential,  v' ,  turbulent  velocity  components.  The  data 
behind  the  rotor,  Figures  II-C70  -  79,  has  two  striking  features. 


Sharp  discrete  tones  at  multiples  of  2970  Hz.  This  frequency  is 
the  blade  passing  frequency  and  thus  the  tones  represent  th° 
passage  of  the  rotor  blade's  pressure  field  and  wake  and  are  not 
associated  with  the  turbulent  velocity  components. 

The  broader  peaks  at  approximately  4,  10,  13,  and  16.5  KHz  are  the 
turbulent  velocity  peaks. 

The  "cleanest"  spectra  occur  at  the  tip  with  an  increasing  amount  of 
"haystacking"  developing  as  the  hub  wall  is  approached.  In  general  the  peaks 
are  higher  at  the  tip  but  the  overall  level  of  turbulence  appears  h.'gher  near 
the  hub. 

The  two  prominent  peaks  are  most  likely  due  to  the  vortices  being  shed  from 
the  pressure  and  suction  surfaces  of  the  blades.  There  appears  to  be  a  decrease 
in  the  amplitude  of  these  two  peaks  as  the  immersion  depth  increases.  This 
would  be  expected  since  the  blade  relative  Mach  number  is  decreasing.  It  was 
shown  in  Section  II-B3  that  the  lower  relative  Mach  numbers  produced  lower 
u'  and  v'  turbulence  levels. 

As  far  as  the  increase  in  the  number  of  "haystacks"  at  the  deeper  immersion 
depths  is  concerned  it  appears  that  the  higher  aerodynamic  loading  in  the 
hub  region  is  playing  an  important  part  in  this  phenomenon.  Table  II-C3  shows  the 
design  AP/q  distribution  for  the  D/V  //4  fan.  It  is  evident  that  Lhe  wide 
band  turbulence  (immersion  #5)  and  the  high  loading  region  correspond.  Since 
these  data  are  derived  from  rms  narrowband  results,  the  peaks  in  the  turbulent 
spectra  represent  the  "average"  of  occurances  on  each  blade.  Thus,  it  can  be 
assumed  that  the  "averaged"  turbulence  spectra  are  characteristic  of  each 
blade.  It  then  appears  that  if  noise  and  u',  v*  turbulence  are  related, 
as  reported  in  Section  II-B3,  that  the  hub  region  of  the  rotor  blade  is  as 
important  in  generating  broadband  noise  as  the  tip  region. 

Figures  II-C80  through  85  show  the  v',  u'  turbulence  spectra  behind  the  OGV. 

For  these  traverses  the  hot  film  probe  was  positioned  about  2  chords  downstream 
and  in  between  two  OGV  wakes.  At  this  distance  from  the  rotor  fewer  and  weaker  har¬ 
monics  of  the  blade  passing  frequency  can  be  observed.  The  two  turbulent  peak.', 
at  4  and  10  KHz  are  still  strong.  In  addition  the  u'  spectra  show  a  high 
turbulence  level  at  the  low  frequency  (less  than  2  KHz)  near  the  outer  wall. 

This  is  due  to  the  tip  wall  boundary  layer  which  is  thicker  behind  the  OGV 
than  behind  the  rotor. 

Flour's  II-C86  and  87  show  the  v'  and  u'  turbulence  spectra  for  a  tip  relative 
Mach  number  of  .729  as  measured  behind  the  rotor  at  immersion  #1.  The  two  peaks 
at  about  4  and  10  KHz  can  still  be  observed  although  they  are  lower  in  amplitude 
than  those  at  the  higher  Mach  number  (Figures  II-C70  and  75). 
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Radius,  in. 

A  P/q 

16.71 

.252 

16.15 

.295 

15.55 

.337 

14.91 

.376 

14.28 

.412 

13.51 

.447 

TABLE  II-C3 

Radial  Design  Point  AP/q  Distribution  For  D/V  #4 
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3-D  ROTOR  VASE  NAAROWBAND  ANALYSIS 
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FIGURE  II-C71.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT 
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FIGURE  II-C72.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  v*  AT  IMMERSION  #3  -  D/V  #4 
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FIGURE  II-C73.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  v*  AT  IMMERSION  #4  -  D/V  #4 
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FIGURE  II-C74.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  v'  AT  IMMERSION  #5  -  D/V  #4 
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FIGURE  II-C75.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  u*  AT  IMMERSION  #1  -  D/V  #4 


FUNDAMENTAL  FREQUENCY 


FREQUENCY  -  KHz 

FIGURE  1I-C76.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  u’  AT  IMMERSION  #2  -  D/V  #4 


FUNDAMENTAL  FREQUENCY  2  ND  HARMONIC 


IMMERSION  #3  -  D/V  #4 


FIGURE  II-C78.  3-D  ROTOR  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  u'  AT  IMMERSION  #4  •  D/V  #4 


FREQUENCY  -  KHz 

FIGURE  II-C80 .  3-D  OGV  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  V'  AT  IMMERSION  #1  -  D/V  #4 


IMMERSION  03  -  D/V  04 


IGURE  XI-C82 .  3-D  OGV  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  v’  AT  IMMERSION  #4-  D/V  #4 


FIGURE  II-C83.  3-D  OGV  WAKE  NARROWBAND  ANAIYSIS  OF  TURBULENT  COMPONENT  u»  AT  IMMERSION  #1  -  D/V  #4 
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FIGURE  II-C84.  3-D  OGV  HAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  u*  AT  IMMERSION  #3-  D/V  #4 
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FIGURE  II-C85 


OGV  WAKE  NARROWBAND  ANALYSIS 


IGURE  II-C87  3-D  OGV  WAKE  NARROWBAND  ANALYSIS  OF  TURBULENT  COMPONENT  u'  AT  IMMERSION  #1  -  D/V  #4 


A  comparison  of  Figures  II-C67  and  II-C68  in  the  two  dimensional  cascade 
II-C86  and  87  shows  a  correspondence  at  4  KHz;  however  the  higher  frequency 
peak  occurs  at  12  KHz  in  the  cascade.  Table  II-C4  contains  a  comparison  of  the 
levels  of  the  high  and  low  frequency  peaks.  In  general  the  level  agreement 
is  within  reason. 

Further  research  will  be  required  to  make  a  direct  connection  between  far- 
field  noise  and  the  hot  film  turbulence  measurements.  The  cascade  results 
reported  in  Section  II-C4;  however,  do  indicate  a  direct  connection  between 
turbulence  and  noise.  More  two  dimensional  and  rotating  wake  data  will  be 
needed  as  input  to  the  broadband  noise  theoretical  models  outlined  in  Section 
II-B.  This  type  of  investigation  should  be  pursued  in  the  future. 


Frequency  Cascade*  D/V  #4** 


u' 

v' 

u’ 

v’ 

4  KHz 

.058 

.030 

.080 

.060 

10  or  12  KHz 

.040 

.046 

.042 

.056 

*  Incidence  angle  8° 

Tip  Relative  Mach  Number  .7 

**  Tip  Relative  Mach  Number  .729 
+  Measured  in  Volts 


TABLE  II-C4 

Comparison  of  Peak  Turbulence  Level  Amplitudes  for  the  Cascade  and  D/V  //4 


NOMENCLATURE  -  SECTION  II-C 


d  Diameter  of  cylinder 

f  Observed  vortex  shedding  frequency 

1  length  of  cylinder 

M  Upstream  flow  Mach  number 

Re  Reynolds  number 

Sh  Strouhal  number 

u'  Fluctuating  velocity  in  the  flow  direction 

u'/U  Turbulence  intensity  in  the  flow  direction 

U  Local  mean  velocity  in  the  flow  direction 

Uq  Undistrubed  (upstream)  mean  velocity 

V'  Traverse  fluctuating  velocity 

W'  Spanwise  fluctuating  velocity 

W'/U  Turbulence  intensity  perpendicular  to  the  flow  direction 

x  Downstream  coordinate  in  the  flow  direction 

y,z  Coordinate  perpendicular  to  the  flow  direction 

c/s  Solidity 

c  Airfoil  Chord 

s  Spacing  between  airfoils 

t  Airfoil  maximum  thickness 

i  Incidence  angle 

oL  Air  inlet  angle 
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(D)  Multiple  Pure  Tones  -  Blade  Passing  Frequency  Sub-Harmonics 


(1)  General  Description  of  the  Phenomena.  One  of  the  most  "intriguing" 
aspects  of  noise  generated  by  trausomc  rotors  is  the  phenomena  Qf  multiple  pure 
tones  (MPT's),  These  strong  pure  tones  have  been  observed  to  be  of  greater 
amplitude  than  the  blade  passing  frequency  (Figure  II-Dl)  at  high  engine  speeds 
and  are,  therefore,  significant  in  terms  of  overall  f an/compressor  noise, 

Several  very  distinct  phenomena  which  characterize  MPT's  are  listed  below: 

,  Detected  only  in  the  inlet  duct 

.  Generated  only  when  the  rotor  tip  Mach  number  is  greater  than  one 

.  Strongest  MPT's  exist  at  frequencies  below  rotor  blade  passing 
frequency  and  occur  at  integral  multiples  of  engine  revolutions 

.  Amplitude  varies  with  engine  speed 

.  Lower  frequency  MPT's  appear  as  engine  speed  increases 

The  first  point  is  exemplified  by  the  probe  narrowband  data  shown  in  Figures 
II-D2  and  II-D3.  The  exhaust  probe  (Figure  II-D2)  shows  only  the  blade  passing 
frequency  tone.  The  inlet  probe  at  the  same  fan  speed  v  f’igure  I1-D3)  shows 
multiple  pure  tones  not  only  below  the  blade  passing  frequency,  but  at  higher 
frequencies  as  well. 

The  second  characteristic  i  as  some  far  reaching  effects.  Figure  I I -DA  shows 
the  occurrence  of  multiple  pure  tones  (MPT's)  in  the  Development  Vehicle  with 
and  without  IGV's  (see  Section  V).  It  is  clear  from  this  pre-  atation  that  the 
MPT's  began  to  propagate  when  the  tip  Mach  number  of  '•he  rotor  became  sonic;  and 
not  when  the  tip  relative  Mach  mmber  exceeds  one.  Since  propagation  of  MPT's 
begins  before  tip  relative  Mach  numbers  become  sonic,  the  occurrence  of  bow  shock 
waves  must  be  incidental  to  MPT  propagation. 

The  third  characteristic  is  ulearly  seen  in  Figures  II-Dl  and  II-D3. 

The  fourth  phenomenon  can  be  seen  in  the  fan  acceleration  tracking  filter 
plot  shown  in  Figure  II-D5,  The  level  of  the  10th  per  rev  tone  (Figure  II-D3) 
is  shown  during  a  slow  accel.  The  rise  in  level  is  clear  1}  shown  with  the  tone 
at  a  high  plateau  above  8,000  rpm. 

Finally,  there  is  a  fifth  point.  Both  Figures  II-D4  and  II-D6  show  that  the 
lower  number  MPT's  appear  as  fan  speed  increa&eo.  Thus,  it  appears  that  the 
propagation  and  generation  of  the  MPT's  is  closely  related  to  the  engine  speed, 

(2)  Some  Theoretical  Considerations.  A  theory  which  examines  MPT's 
will  have  to  conform  to  the  five  phenomena  mentioned  above. 

Consider  the  following: 

In  Section  II-A  an  equation  was  derived  for  the  acoustic  pressure  in  an 
annular  duct, 
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FIGURE  II-D2  EXHAUST  PROBE  HARROW BAND  SPECTRAL  ANALYSIS  INDICATING  THE  ABSENCE  OF  KPT'S 


FIGURE  II-D3  INLET  PROBE  NARROWBAND  SPECTRAL  ANALYSIS  SHOWING  MULTIPLE  PURE  TONES 

8989  RPM  -  50  Hz  B.W. 


FREQUENCY  ,  cps 
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TRACKING  FILTER  NARROWBAND  SPECTRUM  OF  10/REV  SPL  vs.  FAN  SPEED 


FIGURE  II-D6  MULTIPLE  PURE  TONE  0 
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In  equation  (II-D1)  the  expression  for  g  defines  the  well  known  cutoff 
speed.  That  is. 


N  =  f  —  "\/l  -M  2  +  M  ] 
c.o.  In*  z  0  J 


JO  nc 
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(II-D2) 


where  Nc  is  the  cutoff  speed  of  the  fan  in  rounds  per  minute.  For  the 
simple  case  of  zero  flow  (Mz  =  =  0)  equation  (II-D2)  becomes 


*nm  30c 
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Equation  (II-D3)  can  be  rewritten  in  terms  of  fan  tip  Mach  number 
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(II-D4) 


It  has  previously  been  shown  (Section  II-A1)  that  the  eigenvalue,  is  a 
function  of  the  radius  ration  h,  spinning  mode  number  n,  the  radial  node  number 


For  the  case  of  a  source  with  no  stator  the  spinning  lobe  number  is  given 


n  =  nd 


(II-D5) 


Then  eouation  (II-D4  )  becomes 
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In  The  case  of  the  MPT  source  it  may  be  hypothesized  that  there  exists 
in  the  inlet  a  single  wave  form  composed  of  many  harmonics.  Figure  II-D7 
illustrates  this  condition.  Each  of  the  harmonics  would  then  occur  at  a 
frequency  equal  to  a  multiple  of  the  fan  rotational  frequency.  Or  in  other 
words  at  the  per  revs  shown  in  Figures  1I-D1  and  II-D2. 

In  reference  1  the  value  of  (M-j;)c#0>  for  various  modes  and  radius  ratios 
was  presented.  Figure  II-Db  is  a  reproduction  of  Figure  1  of  reference  1. 

It  is  clear  that  for  the  case  being  considered  that  the  fan  tip  Mach  number 
must  exceed  unity  in  order  for  unatten-ated  propagation  of  the  type  repre¬ 
sented  by  equation  (II-D1)  to  take  place. 

This  characteristic  of  the  analysis  satisfies  the  second  point  made  earlier. 
Figure  II-D8  further  supports  the  observations  made  in  Figures  TI-D4  and  II-D6. 
That  is,  that  the  lower  per  revs  (lower  values  of  spinning  lobe  number)  pro¬ 
pagate  unattenuated  at  higher  tip  Mach  numbers. 

Then  equation  (II-D1)  can  be  used  to  predict  the  occurrence  of  MPT’s  with 
respect  to  the  fan  speed  at  which  they  occur  and  their  frequency.  The  ampli¬ 
tude  prediction  can  only  be  accomplished  if  the  amplitude  constant  Aj®,  can  be 
determined. 

Thus  far  evidence  has  shown  that  MPT  propagation  occurs  when  tip  Mach  number 
exceeds  one.  It  appears  that  the  source  of  MPT’s  is  the  rotor's  pressure 
field.  That  is,  the  pressure  pattern  exhibited  by  the  rotor  blades  as  they 
rotate. 

Therefore,  let  it  be  assumed  that  at  a  given  radial  location  the  circumferential 
circulation  distribution  is  made  up  of  arbitrary  rectangular  segments  as  shown 
in  Figure  II-D7. 

Each  segment  constitutes  the  circulation  of  an  individual  blade.  Thus  there 
are  B  such  rectangular  segments  around  the  circumference. 

There  are  three  pertinent  variables  associated  with  the  segments  form: 

.  amplitude 
.  duration 
.  interval 

The  amplitude  is  the  magnitude  of  the  circulation  about  each  blade.  In  general, 
there  is  a  slight  variation  from  blade  to' blade  and  therefore  the  amplitude 
varies  circumferentially.  The  duration  is  the  width  of  the  rectangular  segment 
while  the  interval  is  the  distance  between  segments.  Both  of  these  quantities 
will  vary  with  blade  geometry  and  blade  loading  level. 

The  circulation  distribution  of  Figure  II-D7  represents  a  generalization  of 
the  analysis  of  Section  II-A1.  If  che  duration  of  each  wave  is  shrunk  to  zero 
and  the  intervals  are  all  set  equal  to  2n/B  the  case  of  the  impulse  input  results. 
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FIGURE  II-D7  ARBITRARY  BLADE  LOADING 
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The  wave  form  of  Figure  II-D7  may  be  fxpresned  in  i  reference  frame  re¬ 
lative  to  the  rotor  in  a  complex  Fourier  seriea. 
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If  the  wave  form  of  Figure  II-D7  is  available  at  three  radial  locations, 
then  the  coefficient  cn(r)  may  be  expressed  as  a  second  order  polynomial. 
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The  series  equation  (II-D7)  is  expressed  relative  to  the  rotor.  In  "rder  to 
transfer  it  to  the  absolute  reference  frame  it  must  rotate  with  the  wheel 
speed.  Thus, 
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From  Figure  II-D7  it  ir  "<?er.  that  the  circulation  is 
Therefore,  equation  (II-Db)  can  be  expressed  as 


piecewise  constant. 
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Thus,  equations  (II-D10)  and  (II-D12)  are  equivalent  expressions  for  the 
wave  form  of  Figure  II-D7  in  terms  of  n  harmonics. 

The  circulation  can  be  expressed  as 


f2" 

r(r)  -  Vr  <10  (II-D13) 
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where  V  is  the  net  tangential  velocity  in  the  plane  of  the  rotor.  Therefore 
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Using  equation  (II-D10) 
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Equation  (II-D15)  co'-stitutes  a  boundary  condition  for  the  acoustic  wave 
rrcblem  for  which  equation  (11-01)  is  the  solution. 

Equation  1)  A*  the  acoustic  pressure  may  be  recast  for  the  tangential 

velocity  perturbation  as 
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Since  only  positive  values  of  n  will  exist,  equation  (I1-D15)  can  be  expressed 
as 


Vr  =  i  n  C  (r) 
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for  the  n^  harmonic.  Equation  (II-D17)  can  now  be  equated  to  equation  (II-D36) 
for  the  n,  m  mode 


Cn<-> 


2  R  c  A  K  (A  r) 
o  run  n  nm 


(Tl-Dlfi) 


180 


Proceeding  as  in  Section  II-A1 
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Combining  equations  (II-D1)  and  (II-D19),  the  pressure  is 
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where 
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From  equation  (IT-D21)  the  radial  distribution  of  MPT  amplitude  may  be  cal¬ 
culated  for  each  n  per  rev  if  C  (r)  can  be  determined.  The  term  C  (r)  is  de¬ 
pendent  or>  the  distribution  of  Figure  II-D7.  Unfortunately,  therendoes  not 
exist  an  analytical  method  of  determining  this  pattern.  Further  work  is  needed 
in  this  3rea. 

(3)  Some  Conclusions  on  MPT  Generation.  Relevant  test  data  and  the 
theory  contained  in  Section  II-D2  indicate  that  MPT's  are  generated  and  propa¬ 
gate  when  the  tip  Mach  number  of  the  rotor  exceeds  unitv.  It  is  therefore, 
possible  to  avoid  MPT’s  by  simply  holding  below  one.  This,  however,  exerts 
a  severe  penalty  on  turbofan  and  compressor  designs. 

The  distribution  F'gure  I1-D7  implies  that  the  noise  source  sti-co^th 
associated  with  the  rotor  steady  loading  is  distributed  over  what  anounts  to 
sub-harmonics  of  the  blade  passing  frequency.  If  the  pattern  of  Figure  TI-D7 
is  more  uniform  (the  most  uniform!  '’one  it  ion  being  define*,  as  equal  amplitude 
rectangular  segments  of  constant  width  and  spaced  at  a  2i /B  angular  positions) 
than  more  energy  will  be  fed  Into  the  BPF.  Conversely.,  a  more  random  pattern 
produces  more  MPT  content  and  less  BPF  generation.  Tc  must  be  kept  in  mind. 


however,  that  this  source  of  BPF  noise  is  not  the  only  source;  e.g.,  viscous 
wake  interactions,  and  therefore  changes  in  BPF  amplitude  may  not  be  de¬ 
tectable. 

Figure  II-D8  shows  that  MPT  propagation  may  be  delayed  by  decreasing  radius 
ratio.  This  is  true  only  for  lobe  numbers  below  about  25. 

Adlitional  noise  reduction  procedures  will  have  to  await  a  more  precise 
definition  of  Figure  II-D7. 
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III.  FAN  COMPRESSOR  NOISE  TRANSMISSION  AND  RADIATION  MECHANISMS 


(A)  Sound  Transmission  in  a  Turbomacninery  Duct.  The  idea  that  sound  can 
be  transported  along  a  duct  by  plane  waves  has  been  known  for  some  time. 

The  higher  order  modes  of  transmission  were  at  one  time  regarded  as  a  high  fre¬ 
quency  phenomena  of  secondary  importance.  Tyler  and  Sofrin,  however,  illustrated 
the  importance  of  these  modes  in  the  analysis  of  transmission  and  radiation 
of  sound  in  an  axial  compressor. 

Compressors  and  fan  ducts  of  modern  turbofan  engines  are  potent  sources  of 
discrete  frequency  noise.  Tyler  and  Sofrin  had  shown  that  die  radiated  discrete 
frequency  sound  is  strongly  influenced  by  the  walls  of  the  air  intake  duct, 
and  that  this  effect  is  due  to  the  transmission  of  higher  order  modes  in  the 
duct.  This  type  of  sound  wave  was  called  "spinning  modes".  A  property  which  is  not 
found  in  plane  waves,  but  was  found  for  "spinning  modes"  was  the  cut-off 
property:  each  higher  mode  in  a  uniform  duct  fails  to  propagate  if  the  frequency 
falls  below  a  critical  value,  the  cut-off  frequency. 

The  effect  of  flow  on  higher  order  modes  propagation  were  problems  not 
considered  by  the  early  investigators  of  compressor/fan  transmission  and  rad¬ 
iation.  For  instance,  for  a  sound  wave  emitted  from  a  point  and  traveling  through 
a  moving  fluid  the  sound  wave  is  known  to  experience  changes  in  its  wave  length 
(Doppler  effect).  The  moving  fluid  also  produces  an  asymmetric  pressure  distri¬ 
bution  about  the  source.  In  the  limiting  case  of  a  Mach  number  equal  to  one 
normal  to  the  direction  of  propagation  the  wave  ceases  to  propagate  and  a  zone 
of  silence  results  upstream  of  the  source. 

In  the  case  of  the  spinning  modes  generated  in  fan  and  compressor  operation, 
the  effect  is  somewhat  similar.  That  is,  the  perturbation  pressure  is  distorted 
by  the  axial  and  swirl  velocities  in  the  engine.  In  addition  to  pressure  changes, 
it  can  be  expected  that  the  cut-off  speed  and  decay  rate  below  cut-off  will  be 
altered  by  the  moving  medium. 

Further , directivity  measurements  of  the  noise  radiated  from  a  duct  or  annulus 
to  the  far  field  shows  a  pattern  which  is  a  result  of  reinforcement  and  can¬ 
cellation  of  sound  pressure  as  well. 

In  an  effort  to  better  understand  corapressor/fan  acoustic  wave  transmission 
and  radiation  mechanisms  a  number  of  studies  were  instituted.  In  the  sub-sections 
which  follow  analytical  and  experimental  investigations  will  be  discussed  which 
enable  the  prediction  of  fan/compressor  noise  through  ducts  and  into  the  radiation 
field,  the  sound  transmission  through  blade  rows,  and  the  retraction  eftects  or. 
pure  tones. 

(1)  The  Three-Dimensional  Acoustic  Transmission  and  Radiation  Program.  The 
three-dimensional  acoustics  computer  program  (References  1  and  2)  is  a  unique  and 
versatile  ^ool  to  investigate  spinning  mode  discrete  frequency  sound  .rave  pro¬ 
pagation  in  axirlly  symmetric  cylindrical  or  annular  ducts  of  finite  length  and 
radiation  into  .ree  space,  it  incorporates  a  numerical  solution  of  the  wave 
equation  in  a  cylindrical  coordinate  system  subject  to  boundary  conditions  at  the 
walls  of  the  duct,  e.g.,  a  compressor  or  fan  duct.  A  brief  description  of  the  code 
f ol lows . 


Theory  For  Axially  Symmetric  Geometries  with  Uniform  Medium 
Axial  Flow.  Assuming  small  acoustic  pressures  in  a  homogeneous  medium  and 
ignoring  second  order  effects,  (linear  theory)  the  velocity  potential  in  the 
field  is  described  by  the  following  equation 
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where  c 
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=  acoustic  velocity 
=  time 

=  cylindrical  coordinates  as  shown  in  Figure  III-A1 
=  V  /c  =  steady  Mach  number  of  the  fluid  in  the 
+z  -  direction 

=  steady  velocity  in  the  z-direction 
«  velocity  potential 
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=  -  defines  the  velocity  potential  and  u 

is  the  velocity  norma1  to  the  boundary 


The  sound  pressure,  p,  is  related  to  the  velocity  pountial  by  means  of  the 
expression 
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where  p  is  the  fluid  density,  A  sinusoidal  time  variation  of  the  sound 
pressure,  velocity  and  velocity  potential  is  assumed,  i.e., 
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Then  the  expression  for  the  pressure  becomes 
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If  one  considers  an  acoustic  source  located  at  the  origin  of  the  coordinate 
system,  »he  velocity  potential  distribution  in  the  field  is  as  follows*  (Ref¬ 
erence  3) : 
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source  strength 
2irf  =  the  frequency 
i“  =  -1 

and  the  symbol  *  is  used  to  denote  a  complex  number. 

Relation  (III-A5)  satisfies  the  differential  equation  (1II-A1).  A  series 
of  terms  of  the  same  form  as  equation  (IIT-A5)  is  used  to  obtain  the  solution  when 
boundaries  are  present.  Assume  a  specified  geometry,  as  shown  in  Figure(III-A2) . 
The  continuous  boundaries  are  replaced  by  discrete  source  points  associated  with 
surface  elements,  as  illustrate:!  in  Figure  1II-A3.  Each  source  radiates  in  an 
infinite  medium.  The  pressure  in  the  field  is  the  sum  of  the  contributions  due 
tu  all  sources  which  must  satisfy  specified  boundary  conditions.  The  distri¬ 
bution  of  source  strengths  on  the  surface  is  assumed  to  be  continuous  and 
representable  by  means  of  a  Fourier  series. 

Considering  a  point  £  somewhere  in  the  acoustic  field,  the  velocity  potential 
may  be  expressed  as  fellows: 
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where  S  denotes  the  total  surface  of  sources  and  A(s)  is  the  source  strength 
per  unit  area.  R  and  R*  are  defined  as  follows: 
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Because  the  wave  equation  is  linear,  it  is  satisfied  by  relation  (Ill-All). 
The  objective  is  to  determine  the  unknown  source  strength  distribution  A(s), 
such  that  the  boundary  conditions  are  also  satisfied. 

Boundary  conditions  must  be  defined  for  every  point  on  the  boundary  S. 
These  conditions  may  be  a  statement  of  normal  velocity,  pressure,  or  acoustic 
impedance.  Let  one  assume  that  the  normal  velocity  at  the  boundary  S  is 
given  by 
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where  ro  is  the  number  of  lobes. 

Because  of  the  rotating  pressure  field  shown  in  Figure  (III-A4)  it  is 
logical  to  assume  a  source  of  the  following  type: 

o  .  .I  ...  a'  -iwt  2  i(m0  —tut ) 

Source  strength  per  unit  area  =  A  e  =  A  e 


According  to  the  definition  of  the  velocity  potential,  the  outward  normal  de¬ 
rivative  of  potential  can  be  written  as  follows: 


3$_  (III-A14) 

3N  '  'U 


where 


4  = 


4>  e 


-iut 


2  i(m6-wt) 
$  e 


Substituting  equation  (Ill-All)  into  equation  (111-A14)  and  considering 
point  j  on  the  surface  S  one  gets 
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where  N. 
J 


is  the  unit  vector  normal 


to  the  boundary  at  point  j. 
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FIGURE  III-A4  SIX  -  LOBE  PRESSURE  PATTERN  ROTATING  WITH  ANGULAR  VELOCITY 
RAD. /SEC. 


The  boundary  is  described  by  a  -I  short  rings. Each  ring  id  paired  with 
another  ring  of  uniform  source  strength(with  a  phase  variation  of  The 

boundary  conditions  of  each  ring  will  be  satisfied  only  at  discrete  points  at  the  mid¬ 
point  of  the  ring.  Therefore,  j  in  equation  (II1-A15)  takes  on  a  finite  number 
of  val  cs: 
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In  like  manner,  equation  (111-All)  (for  a  rotating  lobed  pattern)  may  be  exp¬ 
ressed  as  follows: 
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This  can  be  used  with  equation  (III-A4)  to  obtain  the  acoustic  pressure  ac  j: 
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where  p'  -  P.e 
3  3 


i(m0-uit) 


For  the  impedance  condition  at  j,  one  mu!*t  have 


zj  ■  '  ^  V  Vj  •  0 

uj 


(III-A19) 


The  minus  '  is  introduced  here  because  one  is  interested  in  the  ratio  of. pressure 
to  velocity  in  the  fluid,  and  impedance  is  usually  defined  as  the  ratio  tor 
the  lining.  A  positive  resistance  will  then  indicate  energy  absorption  at  the 
surface. 

By  substituting  equations  (III-A16)  and  (III-A38)  into  equation  (.III-A19)  one 
obtains: 
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The  Ajt  values  represent  a  set  of  K  unknowns.  Equation  (III-A16,18,  or  20) 
or  a  combination  of  these  equations,  then  gives  the  reeded  matrix  which  can 
be  solved  for  A^.  Although  equation  (III-A18)  was  written  to  describe  the 
pressure  at  point  j  on  the  boundary,  it  may  also  be  used  to  calculate  the  pressure 
at  any  other  position  in  the  acoustic  field. 

The  equations  described  here  were  programmed  foe  the  General  electric  635 
computer.  The  program  permits  the  calculation  of  acoustic  pressures  at  any  point 
inside  or  outside  the  duct  once  the  boundary  conditions  are  properly  defined.  The 
input  to  the  computer  program  consists,  normally,  of  the  prescription  of  the 
blade  passing  frequency,  the  number  of  spinning  lobes  and  the  detailed  description 
of  the  geometry.  The  boundary  conditions  are  used  in  the  following  manner: 

*  The  pressure  or  the  acoustic  particle  velocity  is 
usually  defined  as  a  function  of  location  in  the 
generation  plane.  Their  values  are  an  input  to  the 
transmission  problem  and  are  to  be  ofccained  from  :he 
noise  generation  analysis  (reference  4). 

*  If  one  deals  with  a  hard  wall  duct  (no  acoustic  treat¬ 
ment),  the  particle  velocity  normal  to  the  wall  is 
prescribed  equal  to  zero. 

The  program  is  flexible.  The  boundary  conditions  can  be  varied  and  could  reflect, 
for  example,  a  gradient  in  generation  pressures  from  hub  to  tip. 

(?.)  Experimental  Verification.  In  order  to  validate  the  analysis,  theory- 
experiment  comparisons  were  made  for  the  test  configuration  of  a  phased  speaker 
arrev  and  lor  a  26-inch  diameter  single-stage  research  compressor. 

(a)  Phased  Speaker  Array.  In  order  to  validate  the  analysis  exper¬ 
imentally,  a  test  configuration  was  run  in  a  phased  speaker  array.  This  device 
consists  of  a  wooden  pipe  with  24  sound  sources  located  on  the  outer  periphery 
of  a  35"  diameter  cylinder  a’.JL&  of  the  device.  The  sources  can  be  phased  in  such 
a  way  as  to  produce  rotating  lobed  patterns  similar  to  those  produced  in  an  actual 
jut  engine.  The  phased  speaker  me  del  ■‘s  shown  In  Figure  III-A5. 

Figure  III-A6  to  9  show  a  comparison  between  the  calculated  and  measured  pressure 
distributions  inside  the  duct  ol'  a  four-iobed  pattern  at  a  frequency  of  700  Hz. 

The  results  show  a  pronounced  standing-wav^,  pattern  because  the  frequency  of  700 
Hz  is  only  slightly  higher  than  the  cutoff  frequency  (630  Hz).  The  calculations 
agree  quite  well  with  the  experimental  results.  Small  differences  near  the 
meut'n  of  the  duct  may  be  caused  by  room  reflections  because  the  model  was  not 
tested  in  an  anechoic  chamber. 
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FIGURE  III-A5  PHASED  SPEAKER  MODEL 
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(b)  26-Inch  Diameter  Research  Compressor.  Various  treated  and 
untreated  duct  configurations  were  experimentally  investigated  on  a  26"  outside 
dxameter,  low  speed  single  stage  research  compressor  having  a  0.85  radius  ratio 
(reference  5) . 

The  research  compressor  with  a  24-inch  .lopg  duct  extension  is  shown  in 
Figure  III-A10.  The  compressor  has  56  blades  and  48  outlet  guide  vanes.  This 
yields  an  8-lobed  spinning  sound  wave  for  the  lowest  possible  circumferential 
mode.  Four  rotational  speeds  are  possible  through  the  electric  motor  drive  and 
high  speed  gear  arrangement  shown  in  Figure  III-A10.  A  schematic  representation  of 
the  compressor  model  used  in  the  calculations  and  the  dimensions  are  shown  in 
Figure  Ill-All. 

Sound  pressure  level  measurements  were  taken  in  the  compressor  duct  and 
in  the  far  field  (0°^  0<  90°), Figure  AI 11-12  depicts  a  comparison  of  calculated 
and  measured  sound  pressure  level  directivities  at  2585  Hz. 

As  one  can  see,  the  analysis  predicts  the  data  fairly  well  for  angles  ex¬ 
ceeding  25°.  It  does  predict,  however,  zero  acoustic  pressure  along  the  engine 
centerline  for  any  spinning  lobe  pattern  and,  therefore,  deviates  from  the 
experimental  data.  Crigier  and  Copeland  (reference  6)  experienced  similar 
difficulties  when  they  tried  to  compare  results  based  on  Tyler  and  Sofrin's 
analysis  (reference  7)  with  experimental  data.  They,  too,  predicted  a  zero 
pressure  along  the  engine  axis.  The  Morfey-Dawson  analysis  (reference  8)  gives 
the  same  results.  It  seems  that  failure  to  predict  finite  acoustic  pressure 
on  the  engine  axis  results  from  the  fact  that  these  different  analyses  do  not 
take  into  account  the  refraction  phenomena  which  arise  from  inlet  and  exhaust 
flow  effects.  In  fact,  refraction  effects  were  found  to  be  dominant  in  establishing 
some  radiated  sound  directivity.  A  special  experimental  investigation  was 
therefore,  set  up  to  investigate  in  detail  the  refraction  phenomena  and  is 
reported  in  Section  III-C. 
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FIGURE  III-A10  26"  DIAMETER  SINGLE  STAGE  RESEARCH  COMPRESSOR 

WITH  THE  ANNULAR  DUCT  EXTENSION 
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Acoustic  source  strength 
Acoustic  velocity 
Base  of  natural  logarithms 
Sound  frequency 

Number  of  circumferential  lobes 

Axial  flow  Mach  number 

Acoustic  pressure 

Instantaneous  acoustic  pressure 

Radial  distance  in  cylindrical  coordinates 

Polar  radius  in  the  cylindrical  coordinates 

Modified  distance  from  source  to  observer  (Equ.  III-A6  and  III-A7) 
Midpoint  of  a  surface  element 
Duct  surface 
Time 

Time-independent,  Particle  velocity  (Equ,  III-A3) 

Normal  particle  velocity  (Equ.  III-A13) 

Time-independent,  Particle  velocity  (Equ.  III-A3) 

Steady  medium  flow  velocity  in  the  z-direction 
Rectangular  coordinates 
axial  distance  in  cylindrical  coordinates 
Acoustic  impedance 
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Axial  angle  in  cylindrical  coordinates 

Wave  number 

Point  in  the  field 

Medium  density 

Velocity  potential 

Instantaneous  velocity  potential 

Angular  frequency 

Angular  velocity 

Subscripts 


Integer  subscripts 
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(B)  Sound  Transmission  Through  Blade  Rows.  Discrete  frequency  components  of 
rotor  noise  and  rotor/stator  interaction  noise  has  beer,  known  for  many  years  to 
propagate  along  an  engine  duct  as  rotating  spinning  modes.  When  the  tra¬ 
veling  waves  encounter  a  blade  row,  i.e.,  inlet  guide  vanes,  rotor  blades, 
outlet  guide  vanes,  etc.,  a  part  of  the  incident  wave  energy  is  transmitted 
through  the  blade  row  and  a  part,  is  reflected.  An  analysis  of  the  sound  transmission 
through  a  blade  row  is,  therefore,  necessary  to  evaluate  which  portion  of  the  sound 
energy  generated  in  a  multistage  compressor  is  radiated  to  the  atmosphere.  An 
analysis  of  the  phenomena  was  undertaken  by  R.  Mani  at  the  General  Electric  Research 
and  Development  Center.  The  analysis  is  based  on  a  Wiener  Hopf  solution  of  the 
wave  equation  for  a  flowing  medium  and  is  described  in  detail  in  Reference  1  A 
program  based  on  Hani's  work  was  written  for  a  high  soeed  computer.  It  permits  the 
investigation  of  sound  transmissions  in  both  stationary  and  rotating  frames  of  ref¬ 
erence  . 

Experimental  Check.  The  sound  transmission  through  stationary  blade 
rows  was  experimentally  investigated  in  a  two-dimensional  cascade  in  the  Laboratory. 

A  sketch  of  the  experimental  set-up  is  shown  in  figure  III-B1.  Air  was  supplied 
through  the  double  plenum  chambers.  A  muffler  section  seventy-two  inches  long,  con¬ 
sisting  of  one-inch  thick  Scottfelt  on  two  eighteen-inch  walls,  was  designed  to 
damp  out  the  upstream  generated  piping  and  valve  noise.  A  sound  source  con¬ 
sisting  of  an  Altec  Lansing  high  intensity  driver  was  located  downstream  of  the 
vanes  and  mounted  on  a  horn  adapter  radiating  through  a  perforated  plate  set 
flush  with  the  duct  wall.  A  Bruel  and  Kjaer  condenser  microphone  was  positioned 
upstream  of  the  vanes.  A  total  of  two  full  vanes  and  two  one-half  vanes  were  used 
in  these  tests.  The  one-half  vanes  were  mounted  on  the  upper  and  lower  duct 
surfaces.  The  vanes  were  spaced  four  inches  between  mean  lines.  A  more  detailed 
description  of  the  test  setup  is  reported  in  reference  2.  The  sound  transmission 
through  the  vanes  was  investigated  at  flow  Mach  numbers,  varying  from  0.0  to  0.9 
at  frequencies  ranging  from  800  to  6300  Hz.  Some  of  the  experimental  data  are 
reported  in  figure  III-B2  with  the  corresponding  analytical  results.  One  can 
see  that  the  agreement  between  experiment  and  theory  is  fairly  good. 

An  analytical  investigation  was  also  undertaken  to  determine  what  parameters 
are  important  in  the  evaluation  of  che  sound  energy  transmission  through  a  blade- 
row.  A  sketch  of  the  analytical  model  is  shown  in  Figure  III-B3  where: 


Oin is  the  angle  between  the  blade  plane  and  the  circumferential 
direction. 

6  is  the  angle  of  incidence  which  can  be  calculated  from  spinning 
mode  theory  (see  references  ,  for  example).  The  incidence  angle 
is  positive  when  measured  in  the  same  sense  as  the  angle  of 
stagger  c.  .  It  is,  therefore,  restricted  to  values  0  <  y  . 
Otherwise,  the  angle  of  incidence  is  negative.  A  positive  Mach 
number  M  <  0  represents  the  case  of  upstream  sound  wave  propagation. 
The  Mach  number  is  negative  (M  <0)  for  downstream  sound  wave 
propagation. 
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FIGURE  III-B3  MODEL  FOR  SOUND  WAVE  TRANSMISSION 

THROUGH  A  BLADE  ROW 
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Effect  of  Cascade  Solidity.  The  solidity, o  =  —  =  airfoil  chord/cir¬ 
cumferential  spacing  was  varied  from  0.5  to  4.0  for  a  given  airfoil  stagger, 
sound  frequency  and  airflow  Mach  number.  The  ratio  of  energy  transmitted 
through  the  cascade  to  the  incident  energy  is  plotted  in  function  of  incidence 
angle  9,  in  Figure  III-E':.  For  different  cascade  solidities,  the  results  indicate 
that  for  values  of  solidity  above  1.0  (which  is  a  lower  limit  from  an  aerodynamic 
viewpoint)  the  solidity  has  only  a  minor  effect  on  the  transmission  of  sound. 

Effect  of  Sound  Frequency.  The  results  of  an  investigation  of  the  fre¬ 
quency  of  sound  on  its  transmission  through  a  stationary  blade  row  are  presented 
in  Figure  I1I-B5.  It  is  shown,  as  one  would  expect,  that  the  large  energy  re¬ 
duction  is  obtained  at  low  frequencies  and  high  upstream  Mach  numbers  (choking 
effect) . 


Effect  of  Wave  Angle  of  Incidence.  The  angle  of  wave  incidence  to 
a  given  cascade  is  a  primary  factor  in  the  evaluation  of  the  sound  transmission. 
Figure  III-B6  shows  some  results  on  the  effect  of  positive  and  negative  in¬ 
cidence  angles,  respectively,  on  the  energy  transmission  through  a  cascade  at 
a  given  frequency.  One  can  see  that  the  energy  reduction  increases  as  the 
magnitude  of  the  wave  angle  of  incidence  is  increased,  whereas  most  of  the  energy 
is  transmitted  at  zero  incidence  angle. 

Effect  of  Airfoil  Stagger.  The  effect  of  airfoil  stagger  was 
investigated  and  some  results  are  shown  in  Figure  III-B7.  It  is  apparent  that 
airfoil  stagger  can  be  a  primary  factor  in  sound  transmission,  particularly  at 
high  upstream  Mach  numbers. 

Effect  of  Flow  Mach  Number.  Figures  1II-B5  to  7  have  shown  that  the 
amount  of  sound  energy  transmitted  is  a  strong  function  of  the  relative  Mach 
number  to  the  airfoil.  One  can  see  that  inlet  guide  vane  choking  can  be  advanta¬ 
geously  used  in  multistage  compressors  as  demonstrated  in  reference  2. 
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FOR  VARIOUS  BLADE  SPACINGS  AND  INCIDENCE  ANGLES 
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FIGURE  III-B5  EFFECT  OF  INCIDENT  SOUND  FREQUENCY  AND  FLOW 

MACH  NUMBER  ON  SOUND  ATTENUATION  BY  A  BLADE 
ROW. 
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ENERGY  REDUCTION,  dB 


DOWNSTREAM 


NOMENCLATURE  SECTION  III-B 


c  Airfoil  chord 

f  Frequency 

M  Flow  Mach  number 

s  Circumferential  Spacing 

am  Angle  between  the  blade  plane  and  circumferential  direction 

(Compliment  of  stagger  angle)  (See  Figure  III-B3) 

0q  The  angle  of  incidence  of  the  incident  wave  (see  Figure  III-B3) 

a  Solidity  (airfoil  chord/circumferential  spacing) . 
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(C)  Refraction  Effect  on  Sound  Directivity.  The  airflow  in  the  ducts  of  the  tur¬ 
bofan  affects  the  effect  speed  of  propagation.  Sound  propagating  ctgainst  the 
airflow  is  effectively  slower,  while  that  propagating  with  the  flov  is  effectively 
faster  than  the  propagation  speed  in  still  air.  So  when  the  noise  generated  by  a 
rotating  blade  row  stage  reaches  the  inlet  or  the  exhaust  plane,  the  propagating 
wave  must  experience  a  change  in  effective  propagation  speed. 

The  transition  occurs  rather  abruptly  in  the  discharge  because  the  airflow 
velocity  gradients  are  contained  entirely  within  the  exhaust  jet  mixing  zone, 
and  this  zone  is  geometrically  restricted  to  a  relatively  small  volume.  In  the 
inlet,  on  the  other  hand,  the  transition  is  distributed  over  a  relatively  large 
volume  more  or  less  according  to  the  distribution  given  by  the  potential  flow 
lines  into  a  sink. 

Therefore,  it  is  to  be  expected  that  refraction  of  the  sound  rays  will  occur 
for  the  noise  propagating  out  th«  exhaust.  Furthermore,  the  amount  of  the  refraction 
will  depend  upon  the  exhaust  velocity  and  the  temperature  of  the  exhaust  stream, 
since  both  affect  the  speed  of  sound  propagation.  Thus,  sound  from  a  high  velocity, 
high  temperature  turbine  exhaust  should  be  refracted  more  than  that  from  a  low 
velocity,  moderate  temperature  fan  exhaust.  Bu' ,  even  in  the  fan  exhaust,  typical 
exhaust  velocities  are  sufficiently  high  for  the  refraction  effect  to  be  a  very 
strong,  if  not  the  dominant  influence  on  the  far-field  directivity  of  aft-radiated 
internally  generated  (turbomachiner\  )  noise. 

For  the  inlet  radiated  noise,  qualitative  evaluation  of  the  potential  flow 
lines  into  the  engine  inlet  indicate  that  there  may  be  some  slight  shifting  of  the 
radiation  pattern  because  of  refraction,  but  the  effect  is  minor  as  compared  with 
the  conditions  that  exist  in  the  exhaust.  The  modal  patterns  that  exist  in  the 
Inlet  duct  appear  to  be  the  dominant  variable. 

The  effects  involved  are  discussed  in  the  following,  based  upon  the  assumption 
that  the  noise  propagates  along  ray  paths  like  light  waves.  The  trends  indicated 
are  that  aft-radiated  noise  is  redirected  by  refraction  toward  the  sideline  ,  so 
that  the  noise  received  by  a  ground  observer  during  an  aircraft  flyover  is  increased. 

(1)  Basic  Mechanism  of  Refraction.  Refraction  of  sound  occurs  between 
two  regions  of  air  if  there  is  a  difference  in  the  velocities  or  a  difference 
in  the  sonic  speed  between  the  two  regions.  This  effect  is  well  known  in  the 
case  of  light  rays,  for  example,  passing  from  air  into  water(see  Figure  III-C1). 

The  index  of  refraction  as  defined  by  Snells  Law  is: 


p  =  iAn~  Q-1- 

sin  Qr 


water 


Cair  =  sPeed  oi~  propagation  in  air 


Cwater=  spetid  propagation  in  water 
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FIGURE  III-C1  SIMPLE  MODEL  OF  SNELL'S  LAW  OF  REFRACTIVE  INDEX  USING 

AIR  AND  WATER  AS  AN  EXAMPLE 


The  same  thing  happens  with  acoustl"  waves  if  the  two  layers  are  both  air 
but  at  different  temperatures,  since  the  speed  of  sound  depends  upon  temperature. 
Likewise,  the  same  thing  appears  with  acoustic  waves  if  the  two  layers  are 
both  air  but  at  different  velocities  since  the  apparent  speed  of  sound  is  the 
sum  of  the  sonic  velocity  within  the  .Tedium  and  of  the  velocity  of  the  medium, 
itself.  (See  Figure  III-C2) .  Morse  and  Ingard  (Reference  1)  have  published 
the  applicable  derivation  for  the  case  of  current  interest,  with  the  results 
as  given  below  (nomenclature  changed): 


sin  0, 


V  =  — - — 

1  sin  0 


+  V„ 


(III-C1) 


If  the  angle  of  incidence,  0^^  ,  is  known,  then  the  angle  of  refraction,  0  ^  , 
can  be  calculated. 

Simplifying  the  equation  for  the  case  in  which  the  air  in  each  region  is  the 
same  temperature  and  the  velocity  in  region  2  is  zero  (ambient),  then: 


sin0  sin0, 

r  i 

where:  M  =  Mach  number  of  the  flowing  airstream. 

It  should  be  noted  that  the  interface  need  not  be  sharply  defined.  If  the 
velocity  or  temperature  gradient  extends  over  an  extended  width,  the  ray  bends 
locally  according  to  the  local  conditions,  so  that  it  may  proceed  on  a  curved  path, 
but  the  exit  refraction  angle  will  be  the  same  as  if  the  interface  were  sharp. 

Also,  of  course,  the  location  at  which  the  ray  exits  will  be  displaced,  depending 
upon  the  width  of  the  transition  region.  At  far  field  distances  of  interest  in 
regard  to  engine  noise,  this  effect  is  of  secondary  importance. 


The  consequences  of  the  above  equations  are  disci.ssed  in  the  following  to 
provide  an  appreciation  of  the  manner  in  which  ref? action  can  be  expected  to 
modify  predictions  made  without  this  effect. 

(a)  For  the  exhaust,  using  the  simpJified  case: 


When 

0  -  0.  - 
l 

then 

0  -  sine 

90 


1 
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and 


c  JL_ 
“  1+M 


0  -  sin  0 


Ir.  this  case,  the  refracted  sound  reaches  the  far  field  in  an  angle  extending 
bet  /eeri  0°  and  0  *  sin"*  (  1  )  (See  Figure  III-C3) . 

r  1+M 

A  shadow  region  occurs  in  the  aft-end,  the  angle  of  the  shadow  depending 
upon  the  Mach  number  of  the  exhaust: 


M 

a 

0 

.25  .5  .75 

1 

01  =90+0 

r  r 

max  max 

= 

180 

143  132  125 

120 

e1 

r 

s 

Angle  measured  from  engine  inlet 

According  to  this,  at  M  =  0.75,  the  shadow  is  shown  in  Figure  (III-C4). 

Thus,  any  sound  ray  leaving  the  engine  between  125°  and  180°  has  been 
refracted  forward  of  125°.  (This  applied  for  the  fan  generated  noise,  not 
necessarily  for  jet  noise.) 

Further  appreciation  of  the  trend  of  the  directivity  pattern  is  jbtained  by 
noting  the  angle  of  refraction  for  various  angles  of  incidence: 


M  =  0.75 


6i 

=  90  +  0  = 

90 

100 

110 

120 
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0  1 
r 

=  90+0 

r 
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99 
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116 
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122 

123 

124 

125 

This  shows,  for  example,  that  all  sound  originally  directed  between  140°  and 
180°  has  been  redirected  to  between  119°  and  125°;  that  sound  energy  within  a 
40°  arc  has  been  collapsed  within  a  6°  arc.  Qualitatively,  this  indicates  that 
the  region  between  110°  and  125°  of  a  polar  noise  plot  should  have  the  highest 
sound  pressure  level. 

The  directivity  index  change  as  a  result  of  this  effect  is  derived  later 
in  this  section. 
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(b)  For  the  inlet,  potential  flow  lines  and  lines  of  constant 
velocity  potential  for  static  test  operation  are  of  the  character  illustrated 
in  Figure  III-C5  (Reference  2) . 

The  dotted  lines  represent  lines  of  constant  velocity  and  the  spacing  between 
them  is  a  measure  of  the  velocity  gradient.  It  is  noted  that  these  lines  are 
non-spherical  so  that  sound  rays  emanating  from  the  inlet  plane  on  the  axis 
will  generally  intersect  them,  at  an  angle  of  incidence  which,  as  sketched 
in  Figure  II1-C6  cause  a  deflection  toward  the  sideline,  away  from  the  axis. 

The  angle  of  incidence,  0^  ,  is  with  respect  to  the  normal  to  the  velocity 
contout  (this  direction  also  being  the  local  direction  of  the  flow  streamline). 
This  can  be  visualized  more  easily  by  comparing  the  process  with  the  refraction 
of  light  from  a  point  source  inside  a  transparent  ellipsoid  such  that  there 
is  an  abrupt  increase  in  the  propagation  speed  at  the  surface.  If  the  ell  psoid 
were  oriented  with  the  major  dimensions  on  the  axis,  refraction  would  tend  to  de¬ 
flect  the  rays  away  from  the  sideline  instead  of  toward  it  as  shown  in  the  above. 
The  contour  lines  of  constant  velocity  are  such,  however,  that  the  trend  is  to 
cause  sidewise  deflection.  And,  of  course,  the  velocity  gradient  is  gradual,  not 
abrupt  change  of  direction  suggested  by  the  sketches  in  Figure  III-C7. 

The  shape  of  the  constant-velocity  contours  changes  rapidly  in  close  proximity 
to  the  inlet  plane.  Thus,  it  is  necessary  to  consider  the  details  of  the  sound 
source  distribution  over  the  surface  of  the  inlet  plane.  This,  combined  with 
the  fact  that  turbulence  in  the  inlet  is  relatively  low  as  compared  with  the 
mixing  zone  in  the  exhaust,  indicates  that  the  far  field  directivity  will  be 
determined  to  a  large  extent  by  the  modal  patterns  in  the  inlet  duct.  Refraction 
effects  appear  to  be  secondary  for  inlet  radiated  noise.  Therefore,  the  remainder 
of  the  discussion  of  refraction  effects  is  concerned  with  aft  radiated  noise. 

(2)  Refraction  Control  of  Exhaust  Radiated  Noise.  Since  the  exhaust 
radiated  noise  is  increased  toward  the  sideline  because  of  the  refraction 
effect  of  airflow-convected  sound  waves,  there  is  the  possibility  that  by 
changing  the  basic  parameters,  some  control  of  this  effect  might  be  achieved 
to  advantage.  Specifically,  if  the  discharge  air  is  shrouded  with  a  gas  with 
higher  effective  speed  of  propagation,  channeling  of  the  sound  by  internal 
reflections  is  possible.  This  effect  has  been  demonstrated  by  E.  Grande  at  the 
University  of  Toronto  (Reference  3)  using  a  cold  nitrogen  jet  to  contain  pure 
tone  noise  injected  into  it  through  a  hypodermic  needle. 

If  the  exhaust  flow  is  surrounded  by  a  shroud  of  higher  velocity  and/or  higher 
temperature  air,  for  example,  there  will  be  a  critical  angle  of  incidence  within 
which  the  sound  will  be  completely  reflected  from  the  interface  at  the  shroud. 

(See  Figure  III-C8). 

Using  equation  (III-C1)  let  : 

Vj  =  Velocity  of  fan  discharge  air 
r\^  -  Sonic  speed  of  fan  discharge  aix 

V£  =  Velocity  of  hot  shroud  air 
C2  -  Sonic  speed  of  hot  shroud  air 
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FIGURE  III- 


One  could  make  V„  >  and  C ^  ,  and  determine  the  critical  angle  of  inci¬ 

dence  ,  eic»  knowing  it  occurs  when  0^  =  90°  : 


sin6 . 

1C 


+  V, 


sin  90e 


+  V„ 


or  sin6 


ic 


c2  +  (v2-v,) 


For  example,  using  temperature  only  to  affect  C2  and  letting  the  two  streams  have 
the  same  velocity  yield: 


So,  for 

T2 

~ 

1 

2 

3 

4 

eic  = 

90 

45.0 

35.2 

U> 

O 

O 

This  shows  that  at  T2/Tj  =  4,  all  sound  with  an  angle  of  incidence  between  30° 
and  90°  will  be  internally  reflected  and  contained  inside  the  shroud  so  as  to 
be  directed  along  the  axis  of  the  jet  within  the  mixing  zone. 

The  angle  of  it.  ience  yielding  complete  internal  reflection  is  shown  in 
Figure  III-C9  as  a  /.unction  of  the  air  shroud  temperature  for  three  cases  of 
interest: 

(a)  =  V_2  (effect  of  temperature  only) 

(b)  =  Cj  and  =  C2  (effect  of  shroud  temperature  when  the  Mach 

number  of  both  streams  is  the  same.) 

(c)  ^2  =  ^2*  ^1  =  0  (maximum  effect,  combining  both  the  temperature 

and  velocity  effects) 
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Case  (b)  above  is  of  most  practical  interest  since  it  is  achievable  by  having 
the  same  exhaust  pressure  ratio  for  both  streams. 

(3)  Predicted  Directivity  Index  Resulting  from  Refraction.  The 
directivity  index  at  angles  is  defined  as: 


where  D  I  =  Directivity  Index,  dB 
5PL  =  Sound  Pressure  Level 
W  =  Sound  Power 

A  =  Area  over  which  the  sound  power,  W, 
is  distributed 

Before  refraction  the  directivity  index  of  the  incident  sound  energy  is: 


D  I..  =  10  log  Wn,  -  10  logA  A.,  -  const 

01  °  0i  °  0i 

After  refraction  it  is: 

D  L  =10  log  W  -  10  logA  A„  -  const 
0r  °  Or  e  Or 

where  A  =  AAQ  is  a  small  increment  of  area  associated 
with  the  direction  0. 
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Then 


D  1  6r  "  D  *61 


+  10  log  ^9r 


W 


-  10  log 


AA6i 


01 


AA 


0r 


”  D  IQi  -  10  log 


AA 

AA 


6i 

6r 


because  W„  =  Wn. 

0r  0i 


(4)  Refraction  from  a  Simple  Engine  Exhaust.  Thus  the  directivity 
index  of  the  refracted  sound  depends  upon  the  directivity  of  the  incident 
sound  rays  and  upon  the  area  into  which  the  incident  sound  has  been  refracted. 
The  incremental  area  of  a  hemispnere  measured  from  the  exhaust  axis,  is: 


2 

AA  =  2tt  R  (sine  0  -  sin6  ,) 
0  Z  1 


So,  the  incremental  area  of  the  hemisphere  for  the  incident  sound  is: 


2 

AA0i  =  2u  R  (sin  6^  ~  slnG  . 


and  for  the  refracted  sound: 


2 

AA  =  2  ttR  (sine  „  -  sin  6  .). 

9r  r2  rl 


Solving  Equation  (III-C1)  for  sin  8  gives  for  the  exhaust  radiated  noise  in 
the  aft  quadrant: 


sin  8 , 


sin  6 


ci  X  V,  -  V 

1  +  1  p 


sin  0 . 
i 


2  33 


DIer  =  DIei  +  20  Log 


C1  V1  ”  Vn  C1 

_J_  - E  sinO.  -  10  log  f 

C  C  Lc 

o  o  J  *- 


(III-C2) 


where 


C  =  Sonic  speed 

=  Exhaust  velocity 

V  =  Plane  velocity 

P 

subscript  1  =  Engine  exhaust 

0  =  Ambient 


234 


The  angle  at  which  this  directivity  index  applies  is  the  refraction  angle, 
er»  given  by: 


8 


r 


sin 


-1 


sin8^ 


sin0 


(III-C3) 


_  .  _  _  0,  ,  then  D  I.  and  0  are  determined.  Note  that  the 

So.  given  D  I  ana  i  ’  6r  r 

angle  as  measured  from  the  engine  inlet, $, is: 


<t>  =  6+90° 


These  equations  simplify  to  the  following  for  the  case  of  fan  discharge  noise 

under  static  test  conditions  (C,  ^  C  ,  V  =0) 

i  o  p 


D  l0r  =  DI0i  +  20  log  (1  +  M  sin  0^ 


0 

r 


^  *  sin8^ 

Sltl  1  +  M  sine. 

l 


The  directivity  index  for  the  simple  case  has  been  plotted  in  Figure  1I1-C10 
for  various  exhaust  Mach  numbers.  If  the  directivity  index  without  flow,  DI0^  , 
can  be  determined  by  analysis,  then  the  directivity  index  with  refraction, 

DI0  ,  is  established. 

In  general,  accumulated  data  on  turbofan  engine  exhausts  suggest  that,  to  a 
first  approximation,  it  may  be  assumed  that: 


DIn,  =  0  for  aft-end  radiated  noise 

Cl 
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So,  then 


DI„  =  20  log  (1  +  M  sime*) 

0r  1 


This  assumption  is  valid  for  a  monopole  source.  Such  an  assumption  is  in 
conflict  with  the  idea  of  lobes  caused  by  the  phasing  relationships  for  the 
modes  generated  by  turbomachinery  blade  rows.  These  phase  relationships  are 
derived,  however,  on  the  assumption  of  wave  propagation  in  a  stationary, 
homogeneous,  non- turbulent  medium.  In  the  actual  case,  the  mixing  zone  of 
the  jet  exhaust  is  very  turbulent,  causing  random  disruption  to  phase  rein¬ 
forcement  and  cancellation  in  the  far  field.  Also,  such  turbulence  introduces 
scattering  of  the  sound  energy  not  only  in  terms  of  geometrical  distribution 
but  also,  in  terms  of  the  frequency  spectral  distribution.  Since  lobes  occur 
as  a  result  of  phase  reinforcement  and  cancellation,  which  is  quite  dependent 
on  the  particular  sound  wave  length,  even  the  frequency  spreading  causea  by 
scattering  has  the  effect  of  causing  the  sound  energy  to  be  of  more  uniform 
spatial  distribution. 

This  is  particularly  true  in  the  sense  of  time-averaged  sound  pressure 
levels.  Instantaneously,  these  lobe  patterns  might  well  exist,  even  after 
refraction, but  they  can  be  expected  to  sweep  back  and  forth  across  a 
far  field  receiver  point  so  that  from  turbofan  engines  with  strong  pure  tone 
content,  the  listener  senses  a  warbling  sound. 

An  exception  to  this  general  rule  may  occur  for  lc  frequency  sound  in  which 
the  wave  length  is  an  appreciable  fraction  of  the  dis Lance  from  source  to 
receiver.  This  is  well  known  from  ground  reflection  cancellation  and  reinforce¬ 
ment  effects.  For  a  hard  surface,  nearly  perfect  reflector,  low  frequency  sound 
pressure  is  doubled  by  reinforcement,  leading  to  6  dB  increase  in  sound  power 
level.  In  this  case,  the  phasing  is  coherent,  and  any  lobe  patterns  as  might 
exist  for  such  low  frequencies  should  be  stable. 

At  high  frequencies,  however,  the  waves  are  scattered  by  air  currents  or 
turbulence  in  a  random  way  so  that  i.he  net  result  is  the  same  as  having  two 
independent,  uncorrelated  sources,  leading  to  only  3  dB  increase  in  sound 
power  level;  such  a  result  can  occur  only  after  periodic  cancellation  and 
reinforcement  of  the  wave  —  the  period  of  this  effect  may  be  very  long 
compared  with  the  period  of  one  acoustic  wave  cycle.  Experience  has  shown 
that,  in  many  cases,  an  averaging  time  of  the  order  of  30  seconds  is  desirable 
to  obtain  repeatability  of  engine  acoustic  data.  This  is  interpreted  as 
an  averaging  of  the  back-and-forth  sweeping  of  the  lobes,  which  tends  to  equalize 
the  far  field  radiation  pattern  insofar  as  aft-end  lobes  are  concerned.  Thus 
the  remaining  effect,  refraction,  becomes  dominant  for  typical  aft-end 
engine  radiated  noise  as  compared  with  any  calculated  initial  directivity  index 
based  on  lobe  patterns. 

(5)  Effect  of  Adding  a  Hot  Shroud  .  When  internal  reflections  occur, 
such  reflected  sound  is  redirected  primarily  along  the  axis  of  the  jet  within 
an  included  angle  essentially  that  given  by  the  jet  exhaust.  That  portion 
of  the  sound  energy  which  is  not  internally  ‘eflected  still  reaches  the  far 
field.  In  passing  through  the  hot  shroud  it  is  refracted  on  entering  and  re¬ 
refracted  on  leaving,  as  indicated  in  Figure  III-C11. 
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The  shroud  thickness  must  be  at  least  one-quarter  wave  length  in  order  for 
ray  analysis  to  apply.  Its  thickness  will  be  so  small,  however,  as  compared 
with  far  field  distances  at  whict  engine  noise  is  measured  that  the  axial 
displacement  inside  the  shroud  can  generally  be  ignored.  Therefore,  the 
Equations  (II1-C2  and  C3)  can  be  used  for  determining  the  effect  of  the  hot 
shroud  on  the  far  field  directivity  where: 


0  -  6.  -  0.  =  t in 

1  1C 


[  C2  +  <VV  ] 


This  has  been  evaluated  '*ith  the  results  as  shown  in  Figure  III-C12  for  a 
case  in  which  the  shroud  is  a.’r  at  various  elevated  temperatures  and  surrounds 
the  fan-discharge  air  which  is  essentially  at  ambient  temperature.  The  dashed 
lines  of  constant  temperature  ratio  in  Figure  III-C12  show  how  the  data  in 
Figure  III-C10  are  modified  by  tue  use  of  the  hot-air  shroud.  For  a  fan 
discharge  Mach  Number  of  1.0,  for  example,  use  of  a  shroud  with  an  absolute 
temperature  three  times  that  of  the  fan  air  temperature  (instead  of  unity) 
yields  a  change  in  the  directivity  index  from  a  maximum  of  6.0  to  a  maximum  of 
3.0,  or  a  net  reduction  of  3.0  dB. 

It  should  be  noted  that  as  the  fan  discharge  Mach  Number  is  reduced  it  is 
necessary  to  use  higher  shroud  air  temperatures  to  obtain  this  beneficial  effect. 
For  example.  Figure  III-C12  shows  that  for  an  exi.aust  Mach  Number  of  0.4,  a  shroud 
temperature  of  twice  the  exhaust  air  temperature  has  no  effect;  this  happens 
because  the  critical  angle  of  incidence  for  internal  reflection  is  so  large 
that  none  of  the  sound  is  reflected.  In  that  case,  however,  further  increase 
of  the  shroud  temperature  will  result  in  a  sideline  noise  reduction. 
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NOMENCLATURES  SECTION  III-C 


c  Speed  of  propagation 

DI  Directivity  Index 

M  Mach  number  of  flowing  airstream  defined  as,  V/C 

SPLa  Sound  pressure  level  at  angle  0 

W  Acoustic  power 

V  Speed  of  flowing  airstream 

0  Wave  angle 

0^  Angle  measured  from  engine  inlet  (refraction,  or  incident  wave  angle) 

^»i 

<f>  Angle  measured  from  engine  inlet 

y  Index  of  refraction 

Subscripts 
i  incident 

ic  critical  angle  of  incidence 

r  refracted 

P  Plane 
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(D)  Propagation  Effects.  In  order  to  interpret  the  results  of  acoustic 
tests  in  the  field,  especially  in  comparing  data  from  one  site  to  another, 
or  ;rom  ground  to  flight  it  is  necessary  to  consider  the  mechanist  ~  associated 
with  the  propagation  of  sound  in  the  atmosphere. 

The  phenomena  considered  in  this  section  are: 

Change  in  amplitude  due  to  inverse  square  law  divergence. 

Absorption  of  energy  by  the  atmosphere. 

Additional  absorption  of  energy  by  the  ground  boundary  layer. 

Reflection  of  sound  waves  by  the  ground. 

All  four  effects  must  be  considered  prior  to  evaluation  of  data  for 
either  the  purposes  of  comparison  and/or  prediction.  The  first  three  will  have 
a  major  effect  on  the  levels  of  sound  energy  reaching  an  observer;  the  fourth 
will  have  significant  effect  on  the  spectral  signature  of  the  source. 

(1)  Inverse  Square  Law.  The  basic  assumption  on  propagation  of 
sound  energy  considers  a  point  source  of  sound  in  free  space  with 
constant  total  power  output.  The  power  flux  (power  per  unit  area)  is  pro¬ 
portional  to  the  amplitude  squared  of  the  acoustic  pressure  wave.  The  total 
power  passing  through  the  surface  of  any  sphere  (with  the  acoustic  source  at  its 
center)  must  be  equal  to  the  total  power  output  of  the  source.  Since  the  surface 
area  at  radius  r  is  proportional  to  t2,  the  amplitude  squared  of  the  acoustic 
wave  is  proportional  to  l/r^  (Inverse  Square  Law) . 

(2)  Atmospheric  Absorption.  The  phenomena  of  atmospheric  absorption 
of  acoustic  energy  are  associated  with  power  lost  due  to  irreversible  work  done 
by  the  wave  on  the  atmosphere.  The  losses  are  due  to  viscosity,  heat  transfer, 
and  molecular  absorption  of  energy.  For  an  acoustic  wave,  these  losses  are 
functions  of  temperature,  humidity,  and  the  frequency  of  the  s^und. 

Reference  1  presents  atmsopheric  absorption  coefficients  for  full  and  1/3 
octave  bands  over  a  wide  range  of  meteorological  conditions,  based  on  a 
combination  of  analytical  models  and  empirical  laboratory  testing,  and  modified 
in  the  high  frequency  range  to  account  for  observed  spectral  shape  as  it 
affects  measurements  taken  with  proportional  bandwidtn  filters. 

These  coefficients  are  still  based  on  the  premise  of  uniform  atmospheric 
properties  along  the  propagation  path,  and  results  will  vary  in  actual  test 
conditions  depending  on  the  extent  of  non-uniformity  present,  or  on  the 
variations  along  the  propagation  path  from  measurements  of  the  controlling 
parameters. 

(3)  Ground  Boundary  Layer  Absorption.  For  sourse  and  receiver,  or 
observer,  both  near  the  ground,  additional  losses  are  observed  which  are  a 
combination  of  effects  of  the  turbulent  boundary  layer  of  the  air  mass  moving 
along  the  ground,  and  the  refraction  of  sound  caused  by  temperature  and  wind 
gradients  along  the  pr  pagation  path.  Reference  2  provides  the  results  of  a 
number  of  empirical  studies  into  these  effects.  For  measurements  or  predictions 
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|  over  relatively  short  distances,  the  acoustic  shadow  effect  shown  in  Figure 

i  1  and  2  can  be  neglected,  but  the  losses  due  to  the  turbulence  of  the  ground 

|  boundary  layer  (illustrated  in  Figure  3,  Reference  2)  are  present  in  the 

real  test  situation  and  must  be  taken  into  account.  For  the  studies  of 
Section  VI -B  of  this  report,  the  levels  of  ground  boundary  layer  losses  shown 
in  Figure  4,  Reference  2  have  been  used. 

(4)  Reflection  of  Sound  Waves  by  the  Ground.  It  is  not  possible  to 
conduct  acoustic  tests  in  the  field  without  a  reflecting  surface  in  the  vicinity. 

Analytical  considerations  of  the  acoustic  field  for  a  point  source  in  the 
vicinity  of  a  reflecting  surface  as  well  documented  (see  for  sample  Reference  3). 

The  reflection  problem  is  ilxustrated  in  Figure  III-D1.  As  indicated  in 
the  Figure  the  problem  may  be  treated  by  considering  the  total  signal  received 
to  be  the  sum  of  the  direct  signals  at  the  receiver  and  at  an  image  receiver. 

For  an  ideal  reflecting  surface  there  is  no  loss  in  amplitude  of  the  reflected 
signal  and  the  reflection  occurs  with  no  phase  shift  in  the  signal.  For  a 
real  surface  (no  diffuse  reflection)  a  certain  portion  of  the  incident  energy 
is  absorbed.  Thus,  the  signal  to  the  image  receiver  has  a  reflection  coefficient, 
Q,  which  is: 

n  =  ^reflected 

x  -  p 

incident 


The  amplitude  of  the  two  components  (direct  and  image)  of  the  signal  will 
be  different  due  to  three  effects: 

Atmospheric  absorption  and  ground  boundary  layer  absorption. 

The  effect  of  square  law  divergence  over  different  path  lengths. 
Absorption  of  part  of  the  incident  waves  by  the  surface. 

The  effect  of  atmospheric  and  ground  boundary  layer  are  negligible  (with 
regards  to  differences  in  amplitude) . 

In  addition,  there  will  be  a  time  delay  between  the  two  signals,  t . 


Ar  _  Ar 
T  c  Xf 

where : 

Ar  =  difference  in  path  length  (reflected  minus  direct) 
c  =  acoustic  velocity 
X  =  wave  length 
f  =  frequency 
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ACOUSTIC  SOURCE 


FIGURE  III-D1  ILLUSTRATION  OF  REFLECTION  PROBLEM 


If  C  is  defined  as  follows: 
r 

C 

r 

where: 

Pq  •  rms  value  of  pressure  at  microphone. 

=  rms  value  of  pressure  due  to  direct  signal 

(i.e.,  signal  that  would  be  received  in  free  field) 

and  with  the  sound  pressure  level  (SPL)  defined  as: 


SPL  -  10  logio  (f^-)  2 

where 

P  =  pressure  (rms  value) 

P  ^  *  reference  pressure 

then  the  value  of  sound  pressrre  level  at  a  microphone  (including  reflection 
becomes : 


SPI.  =  10  log1() 


or 


SPL  =  SPL,  +  10  login  (C  ) 
free  610  r 

field 


i.e.,  10  log^Q  Cr  may  be  interpreted  as  a  delta  to  be  added  to  the  free  field 
sound  pressure  level  in  order  to  account  for  reflection. 

Now,  may  be  written  (Reference  3) 


246 


where  i  denotes  instantaneous  value,  and 


=  rms  value  of  reflected  signal 
=  rms  value  of  direct  signal 

Piri  x  Pdi 

— - - rr-  is  defined  as  the  autocorrelation  function  (R  +  Ar) 

Pr  x  Pd  T  r 

which  is  given  by  Reference  3  as:  * 


(R  +  Ar) 
r»r 


.  fAf 

SinL2f, 


2*  fc  TJ 


2tt  f  T 
c 


where 


A  f  =  bandwidth  of  noise 
f^  =  geometric  center  frequency  of  band 
T  =  time  delay  between  points 

Now  PjVPj  is  (using  the  inverse  square  law  and  reflection  coefficient): 


!i  „  Id  .  Q 

P  11  *  r  2 

d  r 


where  Z  equals  the  ratio  between  the  reflected  path  length  and  the  direct  path 
length. 

Making  the  necessary  substitutions  yields: 

,  ^Q' ,  2Q  sin  'a\  ,  0Ar  x 

Cr  1  *tZ  +  Z  - 71? -  cos  (  —  > 


where 


a  =  2tt  \  2fc  j 

6  *  (If) 
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To  this  point  the  analysis  has  considered  acoustic  waves  travelling 
in  a  homogeneous  medium  (at  rest)  with  reflection  from  a  plane  surface  that 
may  change  the  amplitude  of  the  wave  (Q)  but  does  not  affect  the  phase  of 
the  reflected  wave.  In  general,  the  atmosphere  is  neither  uniform  nor  at 
rest.  In  addition,  the  surface  may  introduce  a  phase  angle  and  it  may  be 
considered  smooth  only  if  X>>  surface  roughness. 

All  of  these  factors  tend  to  create  an  effective  path  length  differential, 
A  Ty  greater  than  the  geometric  difference,  A  r. 

Defining  the  phase  factor  to  be: 


*  « 


Ar 

Ar' 


The  actual  time  delay,  t',  becomes 


or 


Ar 

Xf* 


T  = 


T/4> 


Making  the  substitution  in  the  equation  for  yields: 


ci  Ar 
sin  d>  1 
£  Ar 

(j>  X 


cos  ( 


$  Ar 
<(>  X 


The  condition  for  a  minimum  value  of  C  (null)  is 

r 


,  >  ,  2n  +  1 

Ar  =  X - r — 


n  =  0,  1,  2,  3  ... 


or  for  fixed  geometry: 


4>  - 


Ar 


(2n+l)  X 


null 


2  null  Ar 
(2n+l)  c 


n  =  0,  1,  2. , 
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The  above  relationship  allows  for  experimental  determination  of  the  phase 
factor. 

For  an  ideal  reflector  and  a  homogeneous  atmsophere  <f>  =  1,  and 

f  _  2n  +  1  c 

null  ideal  2  Ar 

Thus,  the  phase  factor  may  be  expressed  as 

^null  actual 

♦  =  f - 

null  ideal 


Summarizing: 


SPL  ,  ,  =  SPL,  ,,  . .  +  10  Log, A  (C  ) 
recorded  free  field  b10  r 


where 


2 

z 


.  ,a  Ars 

sln(» 

V  4>  X  > 


cosl* 


where : 


Q  =  fraction  of  incident  wave  that  is  reflected  (coefficient 
of  reflection) 

2  ’  Vrd 

4r  ’  rc  '  rd 

rr  =  path  length  reflected  signal 
r^  =  path  length  direct  signal 

X  =  wave  length  of  geometric  mean  frequency  of  the  bandwidth  considered 
g  =2ir  /l  +  (  Af/2fc)2 

9  Af 
ct  —  2tt 

c 
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Af  =  bandwidth 


f  =  geometric  mean  frequency 

,  ,  ^null  actual 

and  4>  =  - - 

null  ideal 


Analytical  determination  of  the  phase  factor  (q»)  and  the  coefficient  of 
reflection  is  for  all  practical  purposes  impossible.  The  studies  of  Section 
VI  B  of  this  report  involve  the  determination  of  <j>  and  Q  experimentally  (as 
functions  of  frequency,  surface,  and  geometry). 

For  broadband  noise  analysis  using  proportional  bandwidth  filters 
(such  as  1/3  octave)  the  relationship  Af/2fc  remains  constant.  Thus,  at 
high  frequencies  (small  wave  lengths)  the  last  term  in  Cr  is  damped  out 
due  to  the  1/A  term  in  the  denominator.  As  an  illustration  of  this,  10  log.^ 
C  for  1/3  octave  ■•anda  is  presented  for  an  ideal  reflector  (Q  =  1,<J>  *  1) 
and  for  a  reflecting  curface  with  Q  =  .6  and  $  =  0.9  in  Figure  III-D2. 

For  this  illustration  the  following  geometric  values  were  used:  hs  =  hr  =  12* 
and  L  =  50 ' . 

When  pure  tones  (4  f  +  0)  are  being  analyzed  the  equation  for  C  may  be 
simplified  to: 

cr  ■  1  +  f2  +  2  §  <TT> 


Since  a  -*■  0  and  6  -*■  1. 

In  this  equation,  there  is  no  damping  term  associated  with  the  last  term 
in  Cr.  For  the  same  conditions  used  previously,  10  log  C r  is  plotted 

in  Figure  III-D3  for  pure  tones. 
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figure  hi- 


NOMENCLATURE  SECTION  III-D 

c  Acoustic  velocity 

C  Defined  as,  (P  /P,)^ 

f  Frequency 

f  Geometric  center  frequency  of  band 

P  RMS  value  of  pressure 

P^  RMS  value  of  pressure  due  to  a  direct  signal 

Pq  RMS  value  of  pressure  at  the  microphone 

Q  Reflection  coefficient  defined  as,  P  ,/Pn, 

*  reflected  °incident 

rf  Path  length  of  reflected  signal 

r^  Path  length  ol  direct  signal 

Z  Defined  as,  r  / r , 

r  d 

a  Defined  as,  2tt  (Af/2fc} 

S  Defined  as,  2tt  1  +  (Af /2f ^ 

X  Wave  length 

Af  Bandwidth  of  noise 

Ar  Difference  in  path  length  (r^  -  r^) 

x  Time  delay  between  two  signals 

t"  Defined  as,  t / 

d>  Phase  factor  defined  as,  f  ,  ,/f  ,,  ,, 

’  null  actual  null  idea 
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IV.  EXPERIMENTAL  INVESTIGATION  OF  VANE/ BLADE  NUMBER  EFFECTS 


(A)  Vehicle  and  Test  Descriptions.  An  experimental  investigation  was  pursued 
to  determine  the  effects  of  vane/blade  ratio  on  fan  generated  and  transmitted  noise. 
To  perform  a  meaningful  investigation,  two  aerodynamically  optimized  configur¬ 
ations  were  tested  on  a  fan  vehicle. 

The  vehicle  (Development  Vehicle  I)  used  for  baseline  comparison  purposes 
was  a  scale  model  of  a  modern,  high  tip  speed,  high  bypass  ratio  front  fan  (Fig¬ 
ure  IV-A1) .  The  main-stage  rotor  consisted  of  38  blades,  and  was  35  inches 
in  diameter;  there  were  80  outlet  guide  vanes (OGV's)and  no  inlet  guide  vanes.  The 
inner  termination  of  the  OGV’s  was  a  circumferential  island  which  was  located 
approximately  80%  of  the  span  in  from  the  tutor  tip.  Below  this  inland  were 
the  IGV's,  rotor,  and  OGV's  of  a  small  fan  stage  designed  to  provide  the 
proper  flow  characteristics  along  the  inner  portion  of  the  fan  flowpath.  This 
stage  utilized  less  than  30%  of  the  fan  flow,  and  operated  at  a  very  low  tip 
speed  .generating  very  low  levels  of  noise.  The  main  rotor  design  tip  speed  was 
1388  feet  per  second,  the  design  bypass  ratio  was  6:1  and  the  average  bypass 
pressure  ratio  was  1.617:1. 

A  new  system  of  main-stage  OGV's,  incorporating  only  48  vanes,  was  employed 
on  the  vehicle  (D/V  II)  for  the  second  series  of  tests  (Figures  IV-A2,3,and  4). 

The  chord  length  of  these  OGV's  was  increased  proportionately,  to  maintain 
the  solidity  of  the  original  design;  otherwise,  the  vehicle  was  very  similar 
to  the  baseline  configuration.  Rotor  tip  speed,  pressure  ratio,  etc.,  remained 
as  before.  Figure  IV-A5  compares  the  operating  line  data  actually  measured 
from  the  two  vehicles;  as  can  be  seen,  they  are  almost  identical  in  performance. 

Acoustic  tests  were  carried  out  on  the  two  vehicles  at  a  Peebles,  Ohio, 

Test  Facility.  The  far  field  tests  employed  microphones  arranged  in  two  concen¬ 
tric  arcs,  centered  around  the  fan  axes  (see  Figure  IV-A6).  The  outer  arc, 
of  250  ft.  radius,  consisted  of  15  microphones  placed  in  10  degree  incremental 
angles  from  20  to  160  degrees.  (Angles  are  measured  relative  to  the  Fan  axial 
centerline,  with  the  inlet  equal  0°.)  The  inner  arc(D/v  II  only)  had  8  microphones 
placed  at  a  150  ft.  radius,  at  angles  of  50  through  120  degrees.  Probe 
microphones  were  also  employed  in  the  fan  inlet,  fan  nozzle,  and  fan  exhaust 
(for  readings  at  certain  selected  speeds).  Several  test  runs  were  made  at 
different  fan  speeds.  The  sound  data  obtained  from  the  microphones  was 
recorded  on  magnetic  tape  for  later  reduction,  employing  standard  one-third 
octave  band  and  50  Hz  narrowband  techniques.  The  two  sets  of  reduced  C  *-a 
obtained  on  D/V  I  and  II  are  tabulated  respectively  in  Appendix  A  and  B. 

It  is  noted  that  the  acoustic  comoarisons  between  D/V  I  and  D/V  II  are 
performed  entirely  on  the  basis  of  ,ie  250  ft.  arc  data.  The  150  ft.  arc  data, 
taken  only  with  D/V  II,  was  used  for  substantiation  of  reflection  phenomena  invol¬ 
ving  propagation  of  noise  to  different  distances  from  a  fixed  location  noise  source. 
The  reflection  phenomena  data  is  shown  in  Section  VT-B,  Figures  VI-B67  through 
VI-B70, 
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I rGURE  IV- A2  CONFIGURATION  SKETCH 


IGURE  IV-A3  OGV  FOR  D/V  II  FRONT  VIEW 


FIGURE  IV- A4  OGV  FOR  D/V  II  REAR  VIEW 


100 


(B)  Discussion  of  the  Test  Results.  In  Section  II-A1  an  equation  was 
derived  for  the  acoustic  pressure  of  the  blade  passing  frequency.  For  zero 
steady  duct  flow  that  equation  can  be  expressed  as: 


P 


oo 

E 


n=-°° 


°?  .  _2  in(0-  t)  , .  x  18Z 

I  A  p  C  nw  e  R  (Anmr)e 

ran  n 

m=-°° 


This  equation  has  the  following  mathematical  properties: 

When  3  is  imaginary  the  mode  is  "cutoff" 

Series  covergence  in  m  decreases  the  energy  in  successive  radial  modes 

Series  convergence  in  n  decreases  the  energy  in  higher  order  lobe  numbers 

In  the  case  under  consideration  the  test  vehicles  D/V  I  and  D/V  II  had  38 
blades  and  respectively  80  and  48  vanes.  That  is,  except  for  the  vane  numbers  the 
two  tests  were  done  on  identically  performing  test  configurations.  The  cutoff 
speeds  and  mode  numbers  for  each  configuration  are  shown  in  Table  IV-B1.  Table 
IV-B1  also  shows  that  the  lowest  two  lobe  numbers  generated  in,  respectively,  the 
48  and  80  vane  configurations  are  10  and  38,  and  38  and  42.  From  equation  II-A18, 
D/V  I  should  generate  5  dB  less  noise  at  5750  RPM  and  2  dB  less  noise  at  8500 
RPM.  These  are  within  1  dB  of  the  measured  reductions. 

Figure  IV--B1  shows  the  blade  passing  frequency  (BPF)  power  level  (PWL) 
versus  tip  Mach  number  for  the  two  vane  numbers.  It  is  clear  that  at  tip 
Mach  numbers  less  than  one  that  the  higher  vane/blade  ratio  results  in  less 
noise.  These  data  were  constructed  from  narrowband  analyses  of  far  field 
arc  data  corrected  to  the  source. 

At  higher  tip  Mach  numbers  the  noise  reduction  is  not  as  great  as  at  lower 
speeds.  This  suggests  that  the  cutoff  mechanism  (to  whate  rer  extent  it  exists) 
is  no  longer  functioning  and/or  that  another  noise  source (other  than  rotor  wake- 
stator  interactions  is  making  a  significant  contribution  to  BPF  PWL.  In  either 
case  the  high  vane/blade  ratio  Ls  seen  to  reduce  noise  in  the  important  low 
tip  Mach  number  range  which  is  encountered  at  aircraft  approach  and  power  cutback 
engine  speeds. 

Figure  IV-B2  presents  the  2nd  harmonic  data.  No  major  difference  between 
the  two  configurations  was  expected  (similar  lobe  number  at  the  2nd  harmonic 
frequency).  There  has  been  a  slight  increase  at  low  speeds  and  a  slight  decrease 
at  high  speeds  with  48  vanes. 

The  1/3  octave  power  spectrum  at  low  speed  shown  in  Figure  IV-B3  clearly 
indicates  the  BPF  PWL  decrease  and  that  no  appreciable  broadband  noise  increase 
has  occurred  with  D/V  II. 


Vane  Configuration 

|n|* 

tn 

Corresponding 
Cutoff  Speeds 

48  Vanes 

10 

0 

2320 

10 

1 

3135 

10 

2 

4111 

10 

3 

4648 

48  and  80  Vanes 

38 

0 

8057 

38 

1 

9344 

38 

2 

10313 

38 

3 

11166 

80  Vanes 

42 

0 

8866 

42 

1 

10189 

42 

2 

11183 

42 

3 

12067 

*  | n |  =  38  -  KV  -  Where  V  =  48  and  80 


TABLE  IV- B1 

Fundamental  blade  Passing  Frequency  Cutoff 
Speeds  with  48  and  80  Vane  Stators 


(RPM) 


263 


z  z 
o  o 

H  r-i 

i  ^ 

5  z 
o  o  _ 

M  H  Z 


o  o 
o  o 


5  5  2 

00  O  00 
00  <f 


O  4 


^  is 


anbaua  oNissva  aara 


to 

ri 


H 


w 


(SJLXVM  CI_0I  :9<*)  HP  ‘l*d  DINOWHVH  PQ3 


S 

o 

M 

[i-i 


265 


266 


The  5PF  and  second  harmonic  directivities  or  low  and  highspeeds  are  shown 
in  Figures  IV-B4  through  IV-B7.  The  low  speed  BPF  directivity  (Figure  IV-B4) 
indicates  noise  reduction  at  most  locations  around  the  arc,  particularly  in  the 
aft  quadrant.  The  directivity  of  BPF  at  high  speed  shows,  as  did  the  BPF  PWL 
of  Figures  IV-B1,  little  change  in  BPF  SPL. 

The  second  harmonic  curves,  Figures  Itf-Ba  and  IV-B7,  show  almost  no  change  at 
low  speeds  and  a  slight  increase  in  seconc  harmonic  level  with  80  vanes  at  high 
speed  probably  due  to  the  lobe  pattern (80-2  x  38  =  4)  present  on  this  con¬ 
figuration. 

Figures  1V-B8  through  IV-B11  r±present  narrowband  data  at  5750  rpm  for  110 
and  120  degree  far  field  angles  for  the  two  OGV  configurations.  These  data  have 
not  been  corrected  for  atmospheric  conditions  and  cannot  be  compared  directly. 
However,  the  relation  of  the  tone  content  of  the  spectrum  with  respect  to  the 
broadband  background  noise  gives  an  indication  of  the  tone  reduction  obtained 
with  80  OCV's.  For  example,  on  Figures  IV-B8  and  1V-B10  for  48  vanes  the  BPF 
tone  (at  3600  Hz)has  a  greater  sig  al  (tone)  to  noise (broadband)  ratio  than  for  the 
80  van  vehicle  shown  in  Figures  IV-B9  and  JV-B11. 

In  summary,  the  high  vane  to  blade  ratio  has  been  shown  to  reduce  BPF  noise 
generation,  particularly  at  subsonic  rotor  tip  speeds. 
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FREQUENCY  -  HZ  X  10-3 

NARROWBAND  ANALYSIS  OF  48  OGV  VEHICLE  FARFIELD 


FREQUENCY  -  HZ  X  10"3 

FIGURE  IV-B9  NARROWBAND  ANALYSIS  OF  80  OGV  VEHICLE  FARFIELD  DATA 


VEHICLE  FARFIEIJ)  DATA 


V.  EXPERIMENTAL  INVESTIGATION  OF  ICV  VS.  NO  IGV  FAN 


(A)  Vehicle  and  Test  Descriptions.  Previous  tests  conducted  in  order 
to  determine  the  effects  of  IGV's  on  fan  noise  have  been  done  in  a  rather 
crude  manner.  A  set  of  IGV's  would  be  placed  arbitrarily  on  a  vehicle  designed 
for  operation  without  IGV's, and  acoustic  test  data  would  be  taken,  for  compar¬ 
ison  to  no~IGV  data.  The  fallacy  in  this  procedure  is  that  addition  of  the 
IGV's  alters  the  fan  performance;  the  vehicles  being  compared  thus  differed  in 
performance  as  well  as  IGV  configuration. 

A  new  fan  vehicle  was  designed,  in  order  to  permit  the  determination  of 
the  acoustic  effect  of  IGV's  versus  no  IGV's  on  fans  having  essentially 
the  same  performance.  The  new  fan  (Figures  V-Al  and  2)  was  similar  to  the 
reduced-number-of-OGV's  test  fan,  which  was  used  as  a  baseline.  It  had 
the  same  dimensions,  rotor  blade  and  OGV  numbers,  etc.  Seventy-two  inlet  guide 
vanes  (Figure  V-A3)  were  added.  The  IGV's  extended  to  about  75%  of  the  rotor 
span,  since  the  lower  portion  of  the  fan  rotor  did  not  require  preswirl. 

Also,  since  the  far  greater  part  of  the  noise  is  generated  at,  and  transmitted 
from,  the  upper  portion  of  the  rotor,  the  region  of  the  fan  which  was  most 
important  from  an  acoustic  standpoint  was  affected  by  these  IGV's.  The  small 
inner  stage  again  remained  unchanged. 

In  order  to  achieve  the  same  performance  (and  stall  margin)  as  the  IGV-less 
fan,  it  was  necessary  to  design  a  new  set  of  rotor  blades  (Figure  V-A4)  and 
OGV*s  (Figure  V-A5) ,  and  to  increase  rotor  tip  speed  by  less  than  2.5%.  The  OGV 
was  designed  with  15°  preswirl  at  the  tip  and  the  rotor  redesigned  accordingly. 
Pertinent  comparisons  between  the  two  fans  are  listed  below. 


IGV-less  Fan  (D/V  II)  IGV  Fan  (D/V  III) 


No.  of  rotor  blades 

38 

38 

No.  of  OGV's 

48 

48 

No.  of  IGV’s 

— 

72 

Inner  Stage 

— 

— Common — 

Rotor- IGV  Spacing 

— 

2  OGV  chords 

Des.  Pt.  Rotor  Tip  Mach  No. 

1.275 

1.305 

Des.  Pt.  Tip  Relative  Mach  No. 

1.371 

1.253 

Des.  Pt.  Pressure  Ratio 

1.617 

1.617 

Des.  Pt,  Fan  Weightflow  (lb/sec) 
Diffusion  "D"  factor 

244 

244 

Hub 

.284 

.256 

Pitch 

.505 

.489 

Tip 

Work  Coefficient  - 

.359 

.398 

Hub 

1.849 

1.7663 

Pitch 

0.964 

0.9478 

Tip 

0.6284 

0.6004 
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FIGURE  V-A3  IGV  FOR  D/V  III 
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FIGURE  V-A4  ROTOR  FOR  D/V  III) 


The  key  goal  was  to  maintain  the  blade  loading  as  closely  as  possible.  With 
the  small  increase  in  tip  speed,  the  tip  loading  did  increase,  but  the 
pitchline  and  hub  loadings  actually  decreased  compared  to  the  IGV-less  fan. 

Spanwise,  the  average  hub-tip  loadings  were  essentially  the  same.  With 

regard  to  work  coefficient,  the  higher  tip  speed  of  the  IGV  fan  caused  only  a  slight 

decrease,  as  compared  to  the  IGV-less  fan.  The  net  lo? lings  of  the  two  fans  were  * 

thought  to  be  close  enough  to  provide  a  direct  compa.ison  on  the  basis  of 

equivalent  loading. 

• 

Data  recorded  from  performance  tests  of  the  two  fans  indicates  that  they 
were,  indeed,  almost  identical  In  performance.  Figure  V-A6  compares  the  two 
configurations  on  the  basis  of  fan  weight  flow  versus  percentage  of  fan  speed; 

Figure  V-A7  compares  their  operating  line  data. 

Acoustic  testing,  identical  to  that  described  for  the  IGV-less  vehicle  in 
Section  IV,  was  conducted  on  the  IGV  fan.  The  two  sets  of  data  were  then  com¬ 
pared,  to  determine  the  effects  of  the  IGV's  on  noise  generation  and  transmission. 

The  IGV's  were  installed  in  such  a  marner  that  their  incidence  angle  could 
be  varied.  As  an  extension  of  the  previous  study,  two  further  sets  of  tests 
were  run,  in  which  the  IGV's  were  set  at  "off-design"  angles.  Specifically, 
the  two  settings  employed  were  9  degrees  tc  either  side  of  the  nominal  angle. 

Comparison  of  the  data  from  these  tests  with  that  from  the  nominal  incidence 
angle  tests  allowed  a  determination  of  the  acoustic  effects  of  IGV  incidence 
angle  variations. 

(B)  Discussion  of  Test  Results.  This  test  sequence  had  the  following  objectives: 

To  test  theoretical  prediction. 

To  quantify  noise  generation  changes  with  IGV's. 

To  identify  the  effects  of  opening  and  closing  the  IGV's. 

Io  examine  the  effects  of  IGV's  on  MPT  propagation. 

With  regard  to  the  first  objective,  the  theoretical  work  done  in  Section  II-A 
predicts  the  curves  of  Figure  V-Bl.  The  prediction  shows  the  difference  in 
blade  passing  frequency  noise  generation  for  the  IGV  fan  and  the  no  IGV  fan  at 
various  blade  row  spacings.  At  wide  IGV-Rotor  spacing  (greater  than  about  .75 
IGV  chords)  the  difference  between  the  two  configurations  remains  essentially 
constant  for  both  inlet  nd  exhaust  noise.  The  inlet  noise  shows  less  than  one 
decibel  increase  with  the  IGV  ,  but  a  more  than  two  decibel  increase  is  pre¬ 
dicted  in  the  exhaust. 

Figure  V-B2  shows  the  blade  passing  frequency  (BPF)  power  level  (PWL),  as 
derived  from  1/3  octave  data^versus  tip  Mach  number.  There  is  no  descernible 
trend;  i.e.,  both  configurations  are  producing  tie  same  noise.  Figures  V-B3 
and  V-B4  show  the  corresponding  directivity  of  the  BPF  at,  respectively,  6150 
and  7600  rpm.  At  the  lower  speed  the  aft  quadrant  BPF  Sound  Pressure  Level  has 
increased  with  the  addition  of  IGV's.  The  inlet  level,  in  fact,  shows  a  decrease. 

This  decrease  is  most  likely  due  to  the  attenuating  effect  of  the  noise  passing 
through  the  IGV's.  At  higher  speeds  (Figure  V-B4)  these  effects  have  been  washed 
out  and  both  vehicles  are  producing  approximately  the  same  BPF  levels. 
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CORRECTED  FAN  FLOW  -  W ^Jq/^  (LBS /SEC) 


I 


FIGURE  V-A6  COMPARISON  OF  D/V  II  AND  D/V  III  WEIGHT  FLOW 

vs .  SPEED 


283 


1.7 
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CORRECTED  FAN  FLOW  -  (LBS /SEC) 

FIGURE  V-A7  COMPARISON  OF  D/V  II  AND  D/V  III  OPERATING  LINE  DATA 
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FIGURE  V-B2  COMPARISON  OF  FUNDAMENTAL  PWL  IGV  AND  NO  IGV  VEHICLES 
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Figures  V-B5  :  :d  v-^6  are  the  1/3  octave  source  power  spectra  at  6150  and 
7600  rpm.  Figure  V-B5  shows  the  BPF  PWL  (4  KHz  center  frequency  band)  to  be 
equal;  however  there  is  an  increase  in  broadband  noise  with  IGV's.  This  in¬ 
crease  is  probably  due  to  the  increase  in  turbulence  level  of  the  air  entering  the 
rotor  due  to  the  IGV  wakes.  The  increased  turbulence  results  in  random  fluctuations 
of  rotor  lift  which  manifest  themselves  in  random  (broadband)  noise.  At  higher 
speeds  a  significantly  opposite  effect  has  occurred.  At  high  speed  the  multiple 
pure  tone  (MPT)  content  of  the  tover  spectra  has  decreased  reducing  the  noise  at 
frequencies  below  the  BPF  significantly  (See  Section  II-D) .  This  reduction 
is  due  primarily  to  the  "blocking"  effect  of  the  IGV's.  There  are  also  secondary 
effects  associated  with  the  change  in  rotor  aerodynamic  loading  although  none  of 
these  effects  are  thought  to  be  significant. 

Figures  V-B7  and  V-B3  are  the  resulting  PNL  directivities  for  the  two  speeds. 

At  the  lower  speed  the  aft  arc  PNL  has  increased  while  the  front  is  essentially 
unchanged.  At  high  speed,  however,  (Figure  V-B8)  a  decrease  in  front  quadrant  PNL 
has  occurred  with  no  significant  change  in  the  aft  quadrant.  The  inlet  noise 
decrease  is  important  since  inlet  wall  area  available  for  acoustic  treating  is 
limited  and,  in  some  instances,  an  inlet  acoustic  splitter  may  be  the  only 
method  of  obtaining  significant  inlet  noise  reduction.  In  the  exhaust,  on  the 
other  hand,  there  is  usually  more  treatment  area  available  which  can  minimize 
the  exhaust  noise  increase  at  low  speeds. 

As  previously  noted  (Section  V-A)  the  72  IGV's  were  variable;  i.e.,  their 
stagger  could  be  changed.  Tests  were  run  with  the  IGV's  restaggered  9  degrees 
open  and  9  degrees  closed.  The  closing  of  the  IGV's  decreases  the  rotor's 
aerodynamic  loading;  while  the  loading  increases  as  the  IGV's  are  opened.  Opened  IGV 
(Figure  V-B9)show  small  increases  in  BPF  PWL:  particularly  at  ]ow  speed.  Thus, 
increasing  rotor  loading  has  increased  the  BPF  noise.  The  power  spectra  at 
7900  rpm  shown  in  Figure  V-B10  shows  small  increases  in  broadband  noise  with 
the  IGV's  opened  although  there  appears  to  be  no  MPT  increase. 

Closed  IGV's  had  no  effect  trausonicaily  and  caused  small  decreases  at  low  speeds 
(Figure  V-Bll),  Figures  V-B12  and  V-B13  show  the  1/3  octave  power  spectra  at, 
respectively,  5750  and  6900  rpm.  Within  the  data  scatter  there  is  little  change 
in  noise  at  the  low  speed.  At  the  higher  speed  there  appears  to  be  a  noise 
decrease  over  much  of  the  spectrum.  Figures  V-B14  and  V-B15  are  comparisons 
of  the  9°  closed  IGV  and  the  no  IGV  data.  The  BPF  PWL  shows  no  discernible  increase 
over  the  operating  speed  range.  The  7900  rpm  power  spectra  shown  in  Figure  V-B15 
shows  a  small  BPF  PWL  crease,  a  significant  decrease  in  MPT  level,  and  no 
appreciable  broadband  >  increase. 

In  summary,  the  IGV  configuration  has  been  shown  to  decrease  MPT  levels 
at  the  expense  of  small  increases  BPF  SPL  in  the  aft  quadrant  and  in  broadband 

noise.  Opening  the  IGV's  (loading  the  rotor)  increased  noise  while  closing  the 
IGV's  brought .broadband  noise  levels  down  to  the  no  IGV  levels. 

Test  data  on  D/V  III  is  tabulated  in  appendices  C,  D  and  E  for  the  3  IGV  settings 
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FIGURE  V-B10  COMPARISON  OF  1/3  OCTAVE  POWER  LEVEL  SFECTRA  IGV’s  NOMINAL  AND  IG 

9°  OPEN  -7900  RPM 


IGV's  9°  CLOSED 
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FIGURE  V-B12  COMPARISON  OF  1/3  OCTAVE  POW^rs.  LEVEL  SPECTRA  IGV's  NOMINAL 

9®CLOSED  -  5750RPM 


FIGURE  V-BI3  COMPARISON  OF  1/3  OCTAVE  POWER  LEVEL  SPECTRA  IGV's  NOMINAL  AND  IGV's 

9®  CLOSED  -  6900  RPM 


FREQUENCY  PWL  NO  IGV’sAND  IGV's  9°  CLOSED 


.ISON  OF  1/3  OCTAVE  POWER  LEVEL  SPECTRA  NO  IGV’s  AND  IGV 
iSED  -7900  RFM 


